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ABSTRACT 


Thxs study is directed toward two problems of radiative transfer 
in dusty nebulae; the effects of dust scattering on observable optical 
and infrared parameters and the accuracy of approximate solutions. 

We solve the equation of radiative transfer in a static and homo- 
geneous^ but not necessarily uniform, distribution of gas and dust 
around a central empty core with a point source of energy at its center. 
^Vhile the physics of the gas is well understood, we characterize the 
dust properties by a phenomenological extinction cross section, albedo 
and parameters describing the anisotropy of dust scattering. For 
ultraviolet photons we simultaneously solve ionization equilibrixmi 
equations for the gas, and for infrared photons we calculate a self- 
consistent dust temperature. Ray tracing is used to solve for the 
angular dependence of the intensity. 

The addition of ionized gas adds two serious complications to 
the solution of the equation of radiative transfer. First, the initial 
ionization structure and volume of the Stromgren sphere are known only 
approximately a prion. Second, the gas goes from almost completely 
ionized to neutral very rapidly at the Stromgren radius. Methods of 
handling these two complications while solving the boundary value 
problem of radiative transfer have been developed. 

We find that in most cases variation of the albedo (the ratio of 
the dust scattering cross section to the dust extinction cross section) 
for ultraviolet photons does not have significant effects on optical 
observable parameters, and in all cases these effects are much smaller 
than the effects of the variation of the frequency dependence of the 
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dust absorption coefficient. Increasing the albedo of the dust strongly 
affects the emitted infrared spectrum from the dust by reducing the 
amount of energy the dust can absorb and by reducing the variation of 
the dust temperature in the nebula. When we compare our solutions 
with more approximate solutions of radiative transfer, we find that 
the diffusion approximation, the on-the-spot approximation, and the 
Eddington approximation are very accurate when the assumptions of the 
approximation are met. 
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1 . INTRODUCTION 


Historical Background and Physxcal Ideas 

The study of HI I regions is very old. In the seventeenth 
century Christiaan Heygens described the Orion Nebula. However, the 
emission mechanisms of gaseous nebula were not understood until after 
1925 when quantum theory made it possible to interpret the emission 
line spectra quantitatively. Even then there was some confusion as 
to the nature of forbidden lines. It was not until 192 ? that 
S. I. Bowen correctly identified the [OIII] ^ 959 A doublet 

that had previously been ascribed to the unknown element "nebulium”. 
Since about I960 radio observations have been used to map the free- 
free emission from HII regions. Molecular emission was discovered in 
1963. Since 1965 NASA has sponsored a program of infrared observations 
from high-flying aircraft. In May 1975 "th-® 9 ^ Gerard P. Kuiper 
Airborne Observatory became operational. With this instrument it 
has become possible to make high quality infrared broad band photo- 
metric observations of HII regions. 

HII regions are important to astrophysics for four main reasons: 
First, they are short-lived phenomenon that may yield information on 
star formation. Second, they trace spiral arms and study of their 
properties as a function of location may yield information on the 
formation and maintenance of the spiral arms. Third, they enable the 
determination of elemental abundances and variations of abundances in 
the galaxy. And fourth, there exists the possibility of determining 
the properties of dust grains in HII regions where indepen- 
dent optical and radio observations can be made to determine the 
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structure of the nebula. Points three and four are also valid for 
planetary nebulae. In addition, the study of planethry nebulae yields 
information on post mam sequence evolution of stars and the formation 
of white dwarfs. 

The source of ionizing photons that maintain the ionization in a 

detectable HII region is, in most cases, a hot central star with an 

^ o 

effective temperature T ^ 30^000 K. Photoionization adds energy 
to free electrons that quickly thermalize to an effective temperature 
5000°K ^ T ^ 10,000°K. If the electron density is less than about 

if _3 

10 cm , then these electrons excite low-lying energy levels in the 
ions which depopulate through forbidden lines. Permitted lines in 
the spectra are produced by recombinations , ^Vhen the central source 
IS very hot, high degrees of ionization can be achieved by photo- 
ionization. The ionization structures of the gas are given by a 
balance between photoionization and recombination. At the temperatures 
and densities found in gaseous nebulae, collisional ionizations of 
hydrogen and helium are negligible. 

The infrared radiation comes from two separate sources. First 
IS free-free and free-bound emission from the gas. The luminosity 

of this emission can be predicted from observed Hp or radio luraino- 

2 3 

sities. However, the second source, which is usually lO to 10 
times larger than the predicted infrared emission from the gas, comes 
from the thermal emission of the dust grains embedded in the nebula. 
These grains compete with the gas for the absorption of ionizing 
photons. The energy absorbed by the dust is reradiated as 
infrared continuum radiation at an effective temperature 
given by a steady state balance between absorption and emission. The 
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evidence that the dust grains are mixed with the ionized gas comes 
first from optical photographs that show dark absorption features 
inside of HII regions and second from the infrared maps which are 
usually centered on the radio ( which are measuring free-free emission 
from the ionized gas) maps. 

Internal dust alters the observable parameters of a nebula by 
absorbing and scattering radiation. The dust competes with the gas 


in absorption of ionizing radiation and, hence, alters the ionization 
structures of the gas and reduces the volume of the Stromgren sphere. 

If dust extinction varies with frequency, then the electron temperature 
structure of the nebula will be qualitatively different from nebulae 
with only gas. 

1.2 Observational Data 


HII regions are excited by 0 and B stars or clusters of stars 
of population I . Throughout the nebula hydrogen and helium are singly 
ionized and most of the other elements are primarily in singly or doubly 

3 -3 

ionized states. The gas density varies typically from 10 to 10 cm 

4 _3 

with compact objects having densities as high as 10 cm . The gas 
is highly non-uniform with bright knots, rims and condensations 
observable on most optical photographs. There la also evidence for 
neutral condensations inside of the ionized reiongs. The mass 

2 4 

of ionized gas in HII regions is usually in the range 10 to 10 
where the lower limit is a strong function of observational selection. 
The spectra of HII regions contains HI recombination lines and colli- 
sionally excited forbidden lines of Nil, Oil, OIII and Nelll. 

Planetary nebulae are isolated shells of gas that have been 


ejected recently while their central stars evolved into white dwarfs. 
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^ o 

The central sources are much hotter (T ^ 100,000 k) than HII regxons 

and the gas e:!dxi'bits higher degrees of ionization, including Helll, 

OIIl, Nelll, NelV and NeV. The mass xn the shell is usually in the 

2 h -3, 

range 0,1 to l.OM^ (the densities are in the range 10 to 10 cm ) 
and it IS expanding at 25 km/sec typically. The mean lifetime of 

4 

planetary nebulae is a few times 10 years , 

A common feature that many 'HII regions, planetary nebulae and 
other infrared objects have is the qualitative shape of their infra- 
red spectra. Gin iFigure 1 we have plotted the infrared sjjectra 
of DR21 (a compact HII region), W3(0H)/IIIS8 (an HII 
region), JIGC7027 ( a planetary nebula) , and NGC IO 68 ( a Seyfert 
galaxy) . All of these spectra have a sharply increasing infrared 

flux in the near infrared {x < lOp), a roughly power- law section in 
the mid infrared (lOp. £ X SI lOOp.) and a Rayleigh- Jeans tail in the 
far infrared (x > lOOp) . The similarity of these spectra indicate 
two things. First, a similar mechanism of producing the infrared 
emission is probably operating in all of the objects. Second, the 
power-law portion indicates that for dust emission models, there 
must be a range of effective dust temperature in the model and thus 
that radiative transfer is important in describing the temperature 
variaticaa. The absorption feature at lOji that is apparent in the 
spectra of NGC 7 O 27 and NGCIO 68 is believed to be a feature of 
silicates. Ihis feature is observed in the spectra of most galactic 
infrared sources. 

1.3 The Model 

In order to study and understand the effects of dust extinction 
(absorption and scattering) on the observable optical and infrared 
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S|r(u) (Jy) 



WAVELENGTH 

Fig. 1 . Infrared Spectra of HII Regions. DR 21 is a 
compact HII region, w 3 (OH)IRS 8 is a HII 
region, NGC7O27 is a planetary nebula, and 
NGCI068 IS a Seyfert galaxy. 
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parameters of dusty nebulae and to assess the accuracy of more 
approximate solutions of radiative transfer^ we are going to construct 
a numerical model that not only calculates the ionization structure 
of the gas but also that is able to handle general anisotropic dust 
scattering. We will assume a static and homogeneous (but not 
necessarily uniform) spherical distribution of gas and dust around a 
central empty core with a point source of energy at its center. We 
assume spherical geometry because the radio and infrared maps of many 
nebulae suggest spherical symmetry and also because of the numerical 
simplifications which result from the symmetry. The static and homo- 
geneous assumptions are harder to justify observationally . Dynamical 
models of HII regions (c.f . e.g. Yorke >1977) show that the dynamic 
time scales are much longer than the photon diffusion time scales. 

Our models are, therefore, a "snapshot" in the evolution of an HII 
region. We will ignore the ionization shock fronts which exist at 
the StrSmgren radius (i.e. at the transition zone where the gas goes 
from being almost fully ionized to neutral) . Clumpiness can be taken 
into account in our equations , However, in order to simplify the 
interpretation of our results, we assume the nebula to be homogeneous. 

We will also assume that our model HII regions are just ionization 
bound (i.e. there is no matter outside of the Stromgren sphere) and that 
there is no radiation incident on the nebulae from outside. 

1.4 Method of Solution 

It IS scattering by dust that requires detailed numerical solutions 
of the radiative transfer problem because it makes the equation of 
radiative transfer non-local. We solve the equation of radiative 
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transfer using the quasi-dif fusion method (QDm) described by Leung 
(1975), which is based on methods of solution developed by Auer (I97I) 
and by Hummer and Rybicki (I97I) • 'ni® equation of radiative transfer 
can be solved analytically for stellar photons (i.e. photons from the 
central source which have not interacted with the nebula) because of 
the point source approximation. We therefore find it convenient to 
divide the radiation intensity into two parts. First are stellar 
photons that came directly from the central point source. They remain 
stellar photons imtil they interact with the nebula. Second are 
diffuse photons that are produced in the nebula by scattering or by 
emission. The total intensity is then just the sum of stellar and 
diffuse photons. The QDM is then used to solve for the anisotropic 
diffuse radiation field. 

Leung's method (QDM) first defines the Eddington factor, f (r), 

and a configuration function, ^ (^)j which describe the anisotropy 

of the diffuse radiation field. These functions are used to solve 

a moment or diffusion equation. The functions f (r) and g (r) are 

V V 

initially set equal to the values they would have for isotropic 
radiation. The solution of the moment equation is used to calculate 
the scattering source function and ray tracing is used to solve for 
the angular dependence of the diffuse intensity. With the angular 
dependence, the functions f^(^) ®nd ^y(^) be updated. 

This iterative procedure is continued until convergence of the diffuse 
intensity is achieved. 

The addition of ionized gas adds two serious complications to 
this iterative procedure. First, the initial ionization structure 
and volume of the Stromgren sphere can only be approximated. Thus 
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the iteration procedure will change not only the ionization 
structure but also the outer boundary (or Stromgren radius) of the 
model . The second complication is the very rapid change in the 
ionization at the Stromgren radius. Methods of handling these two 
complications have been a major consideration in the development of 
our numerical models . 

Finally, we will find it convenient to divide the intensity of 

o 

radiation into three frequency bins ultraviolet photons (O < X s 9^2A, 
i.e, photons which can photoionize hydrogen) hereafter refeired to as 
UV photons, optical or OP photons (912A < X < l|i) , and infrared or 
IR photons (Ip i X < °° ) • We separately solve the equation of radia- 
tive transfer in each of these frequency bins, simultaneously solving 
the ionization equilibrium equations for UV photons and calculating 
the dust temperature and infrared emissions for IR photons. 

1.5 Summary of Results 

We summarize here some of the major conclusions of this work. We 
find that the major effect of dust scattering is to transport photons 
( that otherwise would be absorbed in the innei parts of the nebula) 
from the inner to the outer parts of the nebula. Consequently, dust 
scattering can substantially increase the volume of the Stromgren 
sphere over what it would be for a nebula with the same central source 

luminosity and dust extinction coefficient but without scattering. We 
also find that in most cases variation of the albedo (the ratio of the dust 
scattering cross section to the dust extinction cross section) for UV photons 
does not have significant effects on optical observable parameters and in all 
cases, these effects are much smaller than the effects of the variation 
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o± the frequency dependence of the dust absorption coefficient. 

S 

Increasing the albedo of the dust in the frequency ranges primarily 
responsible for dust heating decreases the total infrared luminosity 
of the nebula and the dust temperature at the same relative radius 
by reducing the amount of energy the dust can absorb. Also increasing 
albedo reduces the variation of the dust temperature fiom the inner 
to the outer radius of the nebula and hence makes the emitted IR 
spectrum narrower in frequency. 

When we compare our QDM solutions with more approximate solutions 
of radiative transfer, we find that the diffusion approximation and 
the on-the-spot approximation are very accurate when the assumptions 
of the approximations are met. In addition, the Eddington approxima- 
tion, which assumes that the diffuse radiation field is isotropic, 
gives accurate values for all observable parameters in most cases. 
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2. MATHEMATICAL FORMULATION 


2.1. Equatxon of Radiatxve Transfer 


The equatxon of radxatxve transfer xn spherxcally syimnetrxc nebulae 
xs (c.f. e.g. Chandrasekhar I960) 


fi + 

Br r 

+1 


= - K (r)l (r,ja) 

V V 


( 2 . 1 ) 




+ Js.(r) 


where l^(r,jx) xs the photon xntensxty, 4 xs the cosxne of the angle 
between the photon propagation dxrection and the radxal dxrectxon, 
xs the extinction coefficient, scattering coefficient, 

P^(ju.^ 4 ^) xs the scatterxng phase function which determines the probabi- 
lity of scatterxng from the dxrection into the dxrection p for 

photons of frequency v ^ and J^(^) ^.s the emissxvxty per unit volume, 

time, solid angle, and frequency interval. For UV and OP photons, 
the contribution to J^(^) from recombinations xn the gas and for 
IR photons, 1 (r) xs given by Kxrchhoffs Law as 


Js.Cr) = Sjj^(r,v)NXTj 


V 


( 2 . 2 ) 


where £t-t>(^^''^) infrared emissxvxty of the dust grains and 

IR 


* Throughout this discussion we shall be dealing with number, of photons 
rather than energy of radiation, so that Iy(i*^u) ^.s the number of 
photons per unit area, time, solid angle, and frequency interval. 


10 




(2.3) 


IS the photon emissxon for the dust effective temperature . The 

phase function is normalized so that 


+1 


2j^ = 1 

-1 


(2 A) 


where (Di! is the cosine of the scattering angle. The phase function 
will be expanded in terms of Legendre polynomials as 


1=0 


(2.5) 


where the frequency dependence of the coefficients have been suppressed. 

The normalization condition [Eq, (2.4)] gives od^ = 1 . The average cosine 
of the scattering angle, g^ , and the average of the square of the cosine 
of the scattering angle, g^ , are defined in temis of the a>^*s as 

+1 

( 2 . 6 ) 


+1 

i a>l . gl - I j ’3a a P^(a) 

~1 


+ 1 

15 a)g = gg = da p^(a) 

-1 


(2.7) 
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2.2. Moments of the Intensity 


The n 


th 


moment of the intensity is defined as 

+1 



I • 

V 


( 2 . 8 ) 


Most physical processes in nebulae depend on the mean intensity 
J^(r) = M^(r) . In addition, the flux, H (r) = M^(r) , and the 
K-Integral, = l^(r) are needed to solve equation (2.1). The 

net number of photons passing outward through a spherical shell of 
radius r is 

S^(r) = [3+ir H^(r)] . (2-9) 

Two additional functions will be useful in solving Eq. (2.1). First 
the Eddington factor [l<eung (1975) calls this the anisotropy factor] is 

4 

defined as 


fy(^) = K^(r)/jJ(r) . (2.10) 

This IS also called the variable Eddington factor (Auer and Mihalas, I 97 O) . 
Second is a function introduced by Auer (I97I), 


r 





(2.11) 


where r^ is the inner radius of' the nebula. Leung calls ^ (r) the 
configuration function. The Eddington approximation is equivalent to 
setting fy(^) ~ 1/3 and C^(^) ~ 1 • Throughout this discussion, 
f^(r) and C^(^) will refer to the diffuse radiation field only. 
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The point source approximation for stellar photons is expressed as 

I*(r, 4 ) =2 /(r) , - (2-12-) 

and the moments of the stellar intensity are, 

/(r) = H^(r) = K^(r) . ( 2 . 13 ) 

V V V 

2 . 3 * Combined Moment Equa tion 

The intensity in Eq. (2.1) is first divided into the diffuse and 
stellar parts, noting that only the diffuse radiation field has sources 
in the nebula. Integrating Eq. (2.1) over /j, from -1 to +1 using 
the weighting functions 1 and p. , gives the zeroth and first moments, 
respectively. The zeroth moment equations are 

dS*(r)/dr = - K^(r) S*(r) (2.U) 

for stellar photons, and 


2 dr L V J V v'' V V 


a^(r) 


( 2 . 15 ) 


for diffuse photons, where *^^(i^) and K® (r) are the absorption and 

Si s 

scattering coefficient, respectively (K^ = + K^) . The solution to 


Eq. ( 2 . 1 ^) IS 


S^(.) = a 


T 

- V 


(2.16) 


where S (r^) will be determined from the boundary conditions, and the 
V 0 

extinction optical depth T is defined as 
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r 


T 

V 


(r) 






dr 


r 

0 


(2.17) 


Equation (2.16) combined with Eqs. (2.9), (2.12), and (2.13) completely 
describes the transfer of stellar photons once *^y(^) known. The 
remainder of this chapter will be concerned with solving for the diffuse 
radiation field. 

The first moment equation for diffuse photons is 


h : ? [3 - <(>•)]= 


(2.18) 


Substituting Eqs. (2.10) and (2.H) into Eq. (2.18) and using Eq. (2.I5) 

to eliminate H^(r) gives the combined moment equation 
V 


1 d 

f 

2 

r 

2 dr '* 
r 

C^(r)[Vr)-.j^K=(r)] 





dr 


-g^K;(r)r^ /(r) 
_K^(r)-g^Kj(r) _ 


+ 3 (^) 
\) 


(2.19) 


Subsequently, Eq. (2.I9) will be referred to as the JAY equation. Given 

the source function and the functions f (r) and ^ (^) j this is a 

V V 

D. <, 

linear second-order differential equation for J (r) that is easy to 

V 

solve once the boundary conditions have been specified (see Appendix A). 
Before discussing the rest of the quasi-dif fusion method (QDM), some 
simplifying notation and normalized fluxes will be introduced. (On the firsi 
time through, the reader should skip directly to section 2 . 5 ) 
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2.h. optical Depth, and Normalized Fluxes 


The dust in this model is characterized primarily by i-ts extinction 
optical depth, dT , at the Lyman limit (v = Vq, hy^ = 13.6 eV) as 

dT = dr = dr (2.20) 

where dust number density at the inner surface of the 

nebula, (r) describes the radial variation of the dust density 

(d ,(r ) =1) and C7, ^ is the extinction cross section of the dust 
d' 0 d, 0 

evaluated at the Lyman limit. The total extinction optical depth of the 
dust at the Lyman limit is 


T(r) =n^(r^) dr' (r') , 


0 


which, for a power-law dependence of (2 ,(r) , 




( 2 . 21 ) 


( 2 . 22 ) 


gives 


I^t(HV-l] 

I Rq In (R/Rq) 


a 1 
a = 1 


( 2 . 23 ) 


where 


® - "d<“-0> ‘"d.o ' 


*=0 = 


( 2 . 24 ) 


Once the function d^(r) (or a for power-law density variations) has 
been specified, the dust optical depth is characterized by the parameter 
R^' . For uniform nebulae, the value of R^ has little effect on the 
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pxopertxes of the mbdel as long as (t^ = where is 

the outer radius of the nebula) , However, if CC > 1 , then R^/(a-l) is 
the value of T as r -* co , and the piroperties of the model become 
dependent on the value of R^ , 

The frequency dependence of the dust extinction optical depth is 
defined by the function f^(v) 

One other parameter, the albedo, is needed to completely specify the 
dust properties in this model. The albedo is defined as 


oi = K^/K 
V V V 


(2.26) 


Only the dust number density is allowed to vary with radius so that 
^ ^d 0 constant throughout the nebulae. 

The normalized mean intensity (and similarly, all moments of the 
intensity), denoted by a hat, is defined as* 


I6n J (r) 




(2.27) 


where J has units of 1/Hz and S_ is the number of ionizing photons 
V u 

per second coming from the central point source, i.e. 


■I- 


^0 


Sg = / dv[16^ H* (r^)] . 


(2.28) 


* For nebulae without dust, the factor n^(r)a^j q in Eq. (2.24) and in 
other equations throughout this discussion is^ replaced by n(r)ofQ 
where h(r) is the hydrogen number density and is the hydrogen 

photoionization cross section at the Lyman limit. 
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The normalized net flux through a spherical shell is defined as 


H (r) 



so that for stellar photons. 


u 

/• 


d S (r ) = 1 . 

V V 0' 


V, 


0 


(2.29) 


(2.30) 


In a similar manner, the normalized emission "source fmiction" ( the 

source function is usually defined as j (r)/K (r)) is defined as 

V V 

I6it^ R^ J (r) 

Sj (r) = (2.31) 

%o] 


and the normalized extinction, absorption and scattering coefficients 
are 


K 

V 



cr 


d, 0 


*< = <(-)/”d(-) 

•'v - 


d,0 


d,0 



( 2 . 32 ) 


For nebulae without gas, K = f ^(r) . However, K = 03 m all cases 

' V d' ' ' V V 

because gas does not scatter radiation, resonance scattering of Lyman-cx 
radiation is, in fact, absorption followed immediately by an isotropic 
emission (see §3 9). 
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Although the reason for defining these normalized functions may 


appear obscure at the moment, it will be seen later that once the dust 
parameters are specified, the properties of the model are determined by 
one parameter, namely , The value of this parameter determines the 

total extinction optical depth of the dust and determines the relative 
importance of gas and dust. 

Finally^ in terms of these normalized parameters the JAY equation 
takes the form 


dT 




d 

dT 


' 

[f^(T)C^(T)jJ(T)] ^ 


(2-33) 


a 2 . \ d 

k: r"" j (t) + ^ 

V dT 




St O) 


V 


('T) 


K -g, oi 

V J, V 




where only does not depend on optical depth, or equivalently on 

radius. 


2.5- Ray Equation 

Once the JAY equation has been solved, the QDM uses its solution to 
calculate the scattering source function. This section describes how 
ray tracing is used to solve for the angular dependence of the diffuse 

D / ^ 

radiation field. The diffuse intensity is calculated along a 

ray characterized by its distance of closest approach to the center of the 
nebula (see Fig. 2). This calculation^ is necessary because the diffuse 
intensity is, in general, highly anisotropic. For a similar treatment, 
see Hummer and Rybicki (I 97 I) • 

If we transform into (r,b) coordinates, (where r & b^, then 


18 




Fig. 2. Geometry of Ray Tracing 
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(2.34) 


fi = ± -Jl - b^/r^ 

and due to the symmetry, only values of (r,b) for which z ^ 0 need 

to be considered. Because u is double valued, is separated into 

for 0 ^ u * 1 and I for - 1. ^ u < 0 . I^(l ) denotes the diffuse 
V V ^ s>' 

intensity along a ray going in the outward (inward) direction. By this 
division, only positive values of fi need to be considered. 

Along a ray, Eq. (2,1) for diffuse photons becomes 



nBl'^/Sr = - K (r)!'*' + j (r) + 
V V V V 

S® (r,n) 

(2.35) 

and 

iuai“/Sr = K (r) l” - j (r) - 
V V V V 

2® (r,u) , 

( 2 . 36 ) 


where 2® is the scattering "source function", 
V 


+1 


S" (r,M) 

V 


= h 

2 ^ 


/ 


dti' p^(jLi^M')[l^(r,MO + I_^(i-,ji')] 


( 2 . 37 ) 


and from Pig. 2, it is easy to see that jj, . 

It is now convenient to define an intensity-like function, ^ 

and a flux-like function, V , as 

V 

( 2 . 38 ) 


V 

V 




( 2 . 39 ) 
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Alternately adding and subtracting Eqs. ( 2 . 35 ) ( 2 >^) gives 

jLi &U /Ir = - K (r)v (2.40) 

V V V 

and 


^Svylr = - K^(r)u^ + j_^(r) + 2®(r,|u) 


(2.U1) 


from which V can be eliminated to obtain a second order differential 

V 

equation in U , 

V 




M 


K (r) 
i- V 


Sr 




= 7 ( 2 . 42 ) 


or in terms of normalized parameters 




( 2 . 43 ) 


Subsequently, this shall be referred to as the RAY equation. Like the 
JAY equation, once the source functions are given, this is a linear 
second order boundary value problem that is easily solved for U (r,jLi) 
using the difference equations in Appendix A. From the solution, 

J (r) and all higher moments, f (r) , and C (r) can be calculated. 

The scattering "source function" in normalized units is obtained 
by combining Eq. (2.5) for L = 2 , with Eqs. (2.12) and (2.37) "to give 
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2 




+1 

-1 


= cu. 




(2.44) 




When evaluatxng taken from the solution of the JAY equatic 

and H^(t) and t) are taken from the previous iteration (see § 1.4). 

On the first iteration, K^(t) is given by the Eddington approximation 

V 

j ^D/ \ 

and is set equal to zero. This iteration procedure is explained 

in detail in section 6.2. J^(t) is obtained from the solution of 

Eq. (2.14) combined with Eq. (2.9) giving 


J^(t) = S^(t)/R . 


(2.45) 


Computations have shown that this iteration procedure is rapidly con- 
verging and stable. It will be shown later that the relative strengths 

s E 2 ^ 

of 2^(r,ju) and t (r)/R are determined by the parameter . For 

convenience in calculating fluxes, the RAY equation is solved on an 

(^, h) grid, where experience has shown that the most accurate solutions 

are obtained when the b grid is equal to the r grid. However, ten 

b grid points are added for b < r^ in order to accurately calculate 

moments of I^(r-i(i) at r = r„ . 

V 0 

The moments of I^(r,ju) are calculated from U using the 
relations 
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and 


u- 

=y*dfi j 


jj.) = R ^ / dB TJ^(R,B), (2,46) 


0 


0 


R 


iP(t) = f dll Ijy (R,h) = r"^ / dB B V (R^B) , (2.4t) 

V y V J 


0 


0 


R 


tP(t) = J du 1^ U^(R,^i) = R dB B (R^-B^) U^(R;B) 


0 


(2.48) 


/ \ ''B 

where B = n ,( r„) CT, ^ b . To calculate J and K , it xs assumed 
d'’ 0 d_,0 V V 

that TJ_^^R^BJ varxes Ixtiearly between two adjacent B values, is 

calculated from Eq. (2.40) using a three-point approxxmatxon tor the 
radial derxvative of U_^(R,B) . Again a linear approximation for the 
variation of U^(R, B) with B is used. 


2.6. Boundary Conditions 

In order to make the boundary conditions clear, we will reiterate 
the basic structural assumptions of the model • 

1. Spherically symmetric nebula around a central empty core with a 
point source at the origin. 

2. No radiation is incident from outside of the nebula. 

3. There is no mass outside of the Stromgren sphere. 

The spectrum of the central source must-be specified. For example 

it can be taken to be a black body of arbitrary temperature, or for stars 

o o 

which have an effective temperature in the range 30,000 to 50,000 K, the 
spectra given by the model atmospheres of Auer and Mihalas (1972) have 
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been used, in normalized units, the solution of the equation of radia- 
tive transfer for stellar photons (Eq. is 

-T 

S (t) = S (o) e ( 2 . 49 ) 

V V 

where Eq. (2.30) gives the normalization condition for ®y(0) . This 
form for the stellar flux will allow the Stromgren sphere radius to be 
defined in terms of the normalized flux, , rather than having to 

rely on a less consistent (from model to model) ionization structure 
criterion (e.g. defining the Stromgren radius as the point where hydrogen 
IS 90 ^ neutral). 

The JAY and RAY equations for diffuse radiation are two-point 
boundary value problems with mixed boundary conditions. For the JAY 
equation, the inner boundary condition is expressed physically as the 
fact that there is no source of diffuse photons inside of the central 
core i.e. s"(t = O) = 0 . The less physical boundary conditions on J^ 
are obtained from Eq. (2.18) as 


(t)C (t) jf(r)] 

OT ^ V V S) 


T=0 


= S^(0) 


( 2 . 50 ) 


(jay equation inner B.C.) 


(t)C 


C dT*-' V 

V 


(2. 51) 


T=T, 


^ 1 B 2 >'B, . 


(jay equation outer B.C.) 



where the boundary factor, f^ , 


IS defined as 


• < 2 - 52 ) 

For an isotropic radiation field, f^ - 1/2 . This value is used on the 
first iteration and for the Eddington approximation. For the QDM, 

IS evaluated from the solution of the ray equation. Note should be made 
of the fact that, because U^(r^, /li) = > only the function 

is needed to evaluate f 

V 

The boundary conditions for the RAY equation are physically simple. 

At the symmetry line (r=b) , I = I and at the outer boundary I = 0 . 
Using Eq. (2.40), the inner boundary condition becomes 


and 


where 

IS 


d U /dz ^ = 0 

V r=b 


d U /dzl = 0 

V ' r=r 


0 


d/dz IS used because jLt(r=b) = 0 


for ^ ( 2 . 53 ) 

(ray equation inner B.C.) 

for ^ ^ ( 2 . 5 ^) 

The outer boundary condition 




4 


Bu /3r + K U 


V 


V 


= 0 


r=r. 


(2.55) 


(ray equation outer B.C.) 
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3- PHYSICAL PROCESSES IN THE GAS 


3.1. Basic Assumptions 

> 

All of the physical processes used in this chapter are discussed 

in detail in the excellent book by Donald E. Osterbrock^ Astrophysics 

of Gaseous Nebulae (Osterbrock 197^) hereafter referred to as AGN. We 

review some of the basic processes briefly. For densities found in most 

2 Ij. 2 

gaseous nebulae (i.e. 10 to 10 H atoms per cm'^)^ plasma effects are 
negligible and^ for most atomic processes (except for forbidden transitions' 
electron collisional transitions are unimportant. The gas in this model 
IS composed of n(r) hydrogen atoms per cm"^ , Yn(r) helium atoms (where 
Y m 0.1) and a small amount of trace elements (we include the following 
elements O, N^, C, Ne, Ar and s) . It is assumed^ because the abundance 
of the trace elements is small compared to H and He, that they do not 
affect either the mean intensity or the electron density inside of the 
Stromgen sphere. For computational convenience, it is assumed that all 
photoionization cross sections are zero for frequencies v > 10 . 

Two other basic assumptions are made about the gas. First, all 
atoms are photoionized from the ground state. The mean life-time of a 
hydrogen atom against photoionization, given by 


-1 

T 

ION 


CO 

j dv a^(v) hii 


V 


0 


(3.1) 


where <^„(v) is the photoionization cross section for hydrogen. For 
H 

many nebulae (and all of the model nebulae that we have calculated). 
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3 6 

lO-' sec < ^ sec. However^ the mean life-time for spontaneous 

-k -8 

radxatxve transxtxons of hydrogen excxted states xs 10 sec to 10 sec 

p 

for all states except 2 S for whxch xt is 0.1 sec. All of these mean 

life-times are much smaller than T . The second assumptxon is that 

ION 

all recombinations are to n = 1 or n = 2 states. This xs done for two 

reasons By the assumptions of Case B (AGN § 4 . 2 )^ recombinations to 

states with n ^ 3 eventually lead to population of the n = 2 state. 

Also, the recombination rates to n ^ 3 states are much smaller than the 

rates to n = 1 or n = 2 states. This assumption allows us to reduce 

the number of frequencies that we must consider and thus, to do a more 

detailed and accurate calculation of the important aspects of the problem. 

Finally, the electron temperature can be calculated (see AGN 

Chapter 3 ) and it has been shown (Hjellming I966 and Rubin I968) that the 

variation of T^(^) be quite strong, especially at the Stromgren 

radius. However, the effect of dust is to damp out the variations of 

T (Balick 1975 nnd Sarazin 1977) • this reason, our models assume an 

, 4 o \ 

isothermal electron temperature (T^ =10 K) . This assumption will have 

little effect on the ionization structures of hydrogen and helium because 

- 1/2 

of the weak temperature dependence {f=d ) of the recombination coeffi- 

cients, but will affect the intensities of certain lines. 

3.2. Ionization and Recombination 

The ionization equilibrium equation for any two successive stages 

of ionization 1 and i+l for any element X may be written m terms 

% 

of the total mean intensity, the electron density, n (r), and atomic 

© ^ 

constants as 
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00 


n(X ) 
n(X^) 


; 


d\j cj^i(v) hn j^{^) + 


no: T + 

e ^xH-1 


H^) 


(3.2) 


where 


N 

2 

1=0 


n(x^) = n(x) 


(3.3) 


n(x) IS the total number density of the element X , N is the number 

of ions with photoionization potentials, hv^ , less than lOhV^ , 

cr (v) IS the photoionization cross section of the species X^ , 

X+1 X 

CC IS the recombination coefficient for the reaction X + e -* X + v 

X 

and § and 6 are the charge exchange coefficients for 0 and N 

(AGN p. 36). We have neglected dielectric recombination because, in 

general, it has less than a IC^ effect on the ionization structures 

(Sarazin 1977)^ have neglected collisional ionization because of 

/ n ^ o \ 

the relatively low (T ^10 K « hv /k) electron temperatures found in 
most gaseous nebulae. As was mentioned in §3.1, both 

calculated inside the Stromgen sphere assuming the presence of H , He , 
and dust. Table I gives the photoionization cross sections and the recom- 
bination coefficients for the heavy elements. The frequency integral m 
Eq. (3.2) and all other frequency integrals in this discussion are calcu- 
lated assuming that the integrand varies as a power law between two 
frequency grid points. This assumption is necessary because the ioniza- 
tion edges of the gas result in large variations in J (r) as a function 

V 

of frequency, especially at the Stromgren radius. 
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TABLE I. photoionization CROSS SECTIONst AND 
RECOMBINATION COEFFICIENTS 


TRANSITION 


REF 

^th^^'o 

3 

s 

1 a (xlO ^^cm^/s) 
1 rec^ 

REF 





\mmmmmm 

tJMMIHMH 



0 ° (^p)-o’’ (\) 

2.9*1 

3 

1.002 

2 661 

1 0 



0 ° (^P)-O^ (^D> 

3.85 

3 

1.242 

1.242 

1 5 

3 1 

1 

0 ° (^P)- 0 '*‘(^P) 

2.26 

3 

1.367 

4.311 

1.5 



o'"(^s)-o'^2(^p) 

7.32 

3 

a 583 

3.837 

2 5 

20 0 

1 


3.65 

3 

4 040 


3.0 

51 0 

1 

d*'^(2p)-o‘^'^(^s) 

1.27 

3 

5.693 

0 831 

3.0 

96.0 

1 


0.78 

2 

8.376 

2 6 

3.0 

12 0 

1 

c° (3p)-c'^(®p) 

12.19 

3 

0 828 

3.317 

2.0 

4.7 

1 

C+(2p)_c->-2(ls) 

4.60 

3 

1.793 

1.950 

3.0 

23.0 

1 

c+2(^s)-c'"3(^s) 

1.60 

2 

3.521 

2.6 

3 0 

32.0 

1 


0 68 

2 

4.743 

1.0 

2 0 

75.0 

1 

N (^s)-N‘"(^P) 

n.ii2 

3 

1.C69 

4.287 

2,0 

4.1 

1 

n'(^p)-n‘*'^(^p) 

6.65 ‘ 

3 

2 177 

2.860 

3 0 

22,0 

1 

n’*’^(^p)-'n‘‘'^(^s) 

2.06 

3 

3.489 

1,626 

3 0 

500.0 

1 


1.08 

2 

5.697 

2.6 

3.0 

65.0 


/’"( 2 s)-n‘*' 5 ( 1 s) 

o.M 

2 

7.199 

1.0 

2.0 

150.0 

1 

Nft^(^S)-Ke‘^(^p) 

5-35 

3 

1 586 

3.769 

1.0 

2.2 

1 

Ne*" (^P)-Ne^^(^P) 

4.16 

3 

3-012 

2 717 

1 5 



Ne'*‘(^p)-Ne‘*'^(^D) 

2.71 

3 

3.256 

2 148 

1.5 

15.0 

1 

Ne'*' (^p)-Ne^^(^S) 

0.52 

3 

3.529 

2.126 

1.5 



Ne'^^(^p)-Ne'^^(^s; 

1.80 

3 

4.666 

2.277 

2,0 

44.0 

1 

Ne^^(^p)-Ne’^^(2D 

2-50 

3 

5.060 

2 . 3<^6 

2 5 

91.0 

1 

Ne‘‘^(^p)_Ne‘‘’^(^p' 

1.48 

3 

5.253 

2,225 

2.5 

150.0 

1 
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TABI^ I, (Contd ) 


TRANSITION 

'^^jj(xlO~^®cm^^) 

REF 

''th^'^O 

P 

s 

1 «„,(xlO-^V/s) 

j REF 









N; 3 (^)-/\ 3 p) 

3.11 

3 

7 l 4 l 

1 963 

3.0 

230.0 

1 

Ne’^(^P)-Ng^(^P) 

l.!^0 

3 

9 282 

1.471 

3.0 

280 0 

1 


12.6 

6 

0 762 

1 0 

0 264 



s°(^p)-s'^(^s) 

10.6 

6 

1.471 

1.0 

2.33 

4.1 

1 

s ° (^p)-s'^(^s) 

2,56 

5 

7 . 35 *^ 


3.19 



s’^(’^s)-s’’’^(^p) 

8.20 

6 

1 716 

mm 

0.837 

18 0 

1 

s‘^^(^p)-s'*’^(^p) 

4.37 

6 

2.561 


0.121 

27.0 

1 

s^2(^p)V*^(^s) 

0.290 

6 

3.478 

1 0 

-1.95 

57 0 

1 

s+^(^s)-s‘^ 5 ( 2 gj 

11.8 

6 

5.345 

1.0 

2 24 

120 0 

1 

s'‘‘5(^s)-s'*'^(^s) 

7.56 

6 

6.475 

1 0 

2.23 

170.0 

1 

Ar° (^S)-Ar^(^p) 

30.0 


1.159 

3.0 

2.0 

3.82 

4 

Ar'^(^p)-Ar'*'^(^p) 

30.0 


2.032 

2.5 

2.0 

23.8 

4 

Ar"^^(*'s)-Ar-*-^(^P) 

2.67 


2.995 

2 0 

1.0 

H 8.0 

4 

Ar’^ 3 (^ts)-Ar'^^( 3 p) 

0.591 


4.398 

3 0 

2.0 

84.0 

1 

4 

Ar‘^^(^p)-Ar'*'^(^P) 

0.450 


5.518 

4.61 

2.0 

146.3 

4 

Ar'*'^(^p) -Ar’*'^(^S) 

8.0 


5.693 

0 4 

2.0 

120.0 

4 


t Photoionization Cross Sections 

. »«.IpCv/v^)-= + (1-e) 

(o^j^ and axe the threshold cross section and frequency) 

Beferences 

-1. Aldrovandi and Pe'quignot (1973) 

2. Flower (I 968 ) 

3- Henry (I 970 ) 

4. Kirkpatrick (1972) 

5 . Weisheit ( 197 ^) 

6 . Weisheit and Collins (I 976 ) 
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The ground state recombxnatxon coeffxcxent xs related to the 
photoxonxzation cross section by the equation (agn p. 239 ff) 

hjt ^ . 3/2 


a 


CU O / 1,2 \ 

( 2 # tJ <"^1 


(v) 


h(v^-v)/k 


(3A) 



(v ^ 



where 

cu 

1 

and 

are the statistical weights 

1+1 ) 

and 

M 

e 

and T are the electron mass and 

e 


The frequency dependence of a (v) xs a (v)e 

^x+J. ^x 


-hv/kT 


where cr (v), 

for H and He , is a decreasing function of frequency and hy/kT^ « 1 . 


The spectrum of ionizing radiation from the central source goes approxi- 

2 -hy/kT* * ho 1 , , 

Tnately as v e where T ^ 3 X 10 K . Thus a ® 

X 

and He has a much narrower characteristic width in frequency than does 

1 / N 

the central source and we shall approximate o: 


X X 




( 3 . 5 ) 


where 


‘'eff 


dy y a 


(v)/a^. 


X 


1+1 


1+1 


(3.6) 


and 


a 


,1+1 


03 

■/ 


V. 

1 


X 


(3-7) 


31 



Thxs approxxmatxon xs valxd as long as V - 1 « 1 . We shall 

exi 1 

refer to thxs radxatxon as a ''ixne". 

3-3. Hydrogen 

For hydrogen we defxne the xonxzed fractxon as 


X = n(H'^)/n(H) 

(1-x) = n(H°)/n(H) . (3.8) 
The photoxonxzatxon cross sectxon and recombxnation coeffxcxent (Huiraner 
and Seaton I 963 ) ^re 


'"h(v) = <^0 = ^*3 X 10 


-18 


0 


cm 


^hCv) = 


f(vQ/v) for V 2: Vq 


0 


for V < V 


0 


ttjj = 4.13 X 10 cm^/sec (t^ = 10^ ° k) 


1), o 


(3.9) 


Reconibxnatxon to the ground state xs the only source of UV photons 
from hydrogen, and for T = 10^ °K , equatxons (3.5) - (3.7) gxve 

= 1.06 Vq • The ground states recombxnatxon coeffxcxent xs gxven 
by Eqs. (3.^) and (3-5) wxth cDq = 2 and = 1 . Under the assumptxons 
of § 3 . 1 , there are 3 OP Ixnes produced by recombxnatxons to n ^ 2 states. 
Each recombxnatxon to an excxted state produces a Balmer contxnuum (fia-c) 
line. In a manner similar to the treatment of the n = i recombination 
line, we will treat the spectrum of the Ba-C as a § -function. However, 
because these photons cannot ionize hydrogen, we will set the Ba-C 
effective frequency equal to 0 . 25 Vq - 

Of the recombinations to excited states, two-thxrds eventually 

2 2 
populate the 2 P state and one-third populate the 2 S state. The 
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2 P state depopulates through the emissxon of a Lyman-o; (Ly-cz) photon 
at V = 0-75 Vq • model includes resonance scattering of Ly-a photons 

(.see §3.-9). We assume that this does not al'ter^ however, the ground 
state criterion for photo ionizations. For hydrogen, the resonance 
scattering cross section (at line center) is O'j^(Ly-a) 9300 cr^ for 


it o 2 P 


0 


T = 10 K. The 2 S - 1 S forbidden transition produces a 2 -photon 
e 

continuum with the most probable distribution ~ ~ ''^Ly-CZ * 

The emission rates of these OP photons are 

Sa-C = ^ 


^Ly-a ^Ba-C 

^2-photon = 2 X I 


( 3 . 10 ) 


g 

where a is defined as 

X 


a 


X 


, 1+1 




n=2 


(3.11) 


and cc^ . IS the recombination coefficient to the state n . 


3*^* Singly Ionized Helium 

•J* 

The fraction of He atoms that are singly ionized, y , and the 
fraction that are neutral, y*^ , are defined as 

y*^ =: n(He^)/n(He) = n(He*^)/Yn(r ) 

y^ = n(He^)/n(He) = n(He’’’) /Yn( r) . (3.12) 

In the ionization equilibrium equation [Eq. (3.2)], the ionization poten- 
tial of He*^ IS 2^.6 eV or = I .807 Vq , and the photoionization 
cross section and recombination coefficient (Hummer and Seaton 196 ^) are 
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0 


He 


0 


(v) 




f o('>) = 


He 


-O.T3(v/\Jq - 1,8oT) 


V £ 1.807 Vq 

V < 1 . 8 ot Vq 



0 


T.35 X 10 


-18 


cm 


(3.13) 


a , - 4.31 X 10 cm^/sec (t =10 k) . 
TT + ' e ' 

He 


The addition of helium 10^ hy number) to a nebula has a strong 
effect on the hydrogen ionization structure because H and He compete 
for UV photons with frequencies v ^,1.807 '^q • This effect is lessened 
if T^ ^ 4 X 10 K or if dust preferentially absorbes photons with 
V > 1,807 Vq so that the helium Stromgren sphere is smaller than the 
hydrogen sphere. 

Recombinations of singly ionized helium can produce one of four UV 

lines and one of three OP lines. Ground state recombinations produce 

photons at an effective frequency v „. = I .87 (t = 10 ^ ° K) at a rate 

eii 0 e ‘ 

given by Eqs. (3.4) and (3.5) with 1 and = 2 . Three-fourths 

of the recombinations to excited levels are to the triplet state 2 ^S 

11 q 

and one-fourth are to the singlet states 2 P and 2 S , The 2'^S state 

4 1 

IS metastable with a half-life of i«!=i 10 sec. It decays to the 1 S state 

via emission of a single photon forbidden line at I 9.8 eV (v = 1.456 v^) . 

3 -3 

If ^ 5 X 10 cm y the dominant mode of depopulation is through 

1 1 

collisional excitation to the 2 S or 2 P states. For many HII regions, 

3 -3 

n^ ^ 10 cm so this process will not be considered. Although the 
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3 

lxfe~txme of the 2 S state only margxnally meets the Ixfe-txme criterion 
of § 3 . 1 , a detaxled analysxs of excxted state populatxons xs beyond the 

scope of this work. Of the singlet state recorabxnatxons^ one-thxrd are 

1 11 

to the 2 S state and two-thirds are to the 2 P state. The 2 S state 

depopulates (life-time 0.02 sec) by a two photon continuum with 

hv^ + hv^ = 20. Y eV. The probability that a UV photon is produced is 

0.56 and the most probable frequency of the photon is v = l*l 83 Vq • 

Finally, the 2^P state etaits a He Ly-o: photon at y = 1.522 v^, that is 

reasonantly scattered by helium (at the line center, cr ^-.(Ly-a) ^ 

I, He^ 

^ 4 O , 

2600 a „ for T =10 K) . Each of the excited state recombinations 
discussed above produces a Ba-C line with slightly different frequencies. 

These OP lines are at frequencies v = 0. 2^4-1 , O .285 and 0-351 

4 11 3 

corresponding to recombination to the states 2 P , 2 S , and 2 

respectively. 


3 . 5 * Doubly Ionized Helium 


4“l“ 

The atomic structure of He is hydrogenic with Z 

photoiomzation cross section is, therefore, 

a (v) = CT f (y) 

He^ He , 0 He^ 


= 2 


He’*’ 


(v) ^ 


(4vq/v)- 

0 


V 


< 4v 


0 


0 


r = a A = 1.5T5 X 10 cm^ 
He ,0 ^ 


and the recombination coefficient (Seaton i 960 ) is 

= 2,176 X 10 ^ cra^/sec . 


a. 


He 


++ 


The 


(3.14) 


(3.15) 
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IS defined as 


I * ] ■ 

The fraction of helium that is doubly ionized^ y , 

y'*"*' = n(He^^)/n(He) = n( He"^^ ) /Yn( r) . (3- 16) 

The ionization potential of He^ is 5^-^ eV or ~ ^ Vq • Even the 

hottest 0 stars produce few photons with v ^ ^I-Vq so that doubly ionized 

helium will only play a significant role in the properties of gaseous 

nebulae if T* ^ 10^ " *K (i.e. planetary nebulae). 

Recombinations of He^^ produce four UV lines and no OP lines 

with the assumptions of §3.1. Ground state recombinations produce 

photons at an effective frequency = i4-.06 Vq (T^ = 10^ K) at a rate 

given by Eqs. (3.^) a-ud (3.5) with = 2 and = 1 Reconbinations 

to excited states produce a Ba-C photon with v = . Two-thirds of 

these recombinations result in Ly-a photons with frequency v = 3 vq 

that are resonance scattered by He (cr ^(Ly-a) sa 9300 o ^ at line 

He"^ He’*', 0 

center). One-third of the recombinations result in population of the 
2 

2 S state that depopulates through a two-photon continuum. The most 
probable distribution in frequency is 

He , Ly-a 

emission 'coefficients are given in Table II. 

3.6. Recombination Source Term and Extinction Coefficient 

The recombination source term for UV and OP photons that appears 
in the JAY and RAY equations has the general form 

j (r) = n(r) n (r) fi(r)Ajr ( 3 - 17 ) 

V '' 

where 

n^(r) = l^x + Y(y’*’ + 2y'*"*')j n(r) , (3.18) 
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TABLE II. RECOMBINATION LINES OF HYDROGEN AND HELIUM 


SOURCE 

LINE 

v/v„ 

EMISSIONtt 


l^S Rect 

1.06 

1 


Ly-o: 

0.75 

2/3 

rO 

2^S - l^S 

0.375 

0 1 B 

2 X 3 X «H 

0 

H 

Ba-C 

0.25 

B 

®^H 

« 0 
He 

l^S Rec 

1.87 

«He+ 

He^ 

Ly-a 

1.522 

1/4 X 2/3 

He® 

2^3 - l^S 

1.456 

3A 

He® 

1 

2 S - 1 S 

1.183 

1/4 X V 3 (0 56) 

He® 

3 

2 “^S Rec 

0.351 


He® 

2 ^S Rec 

0.285 

1/4 X 1/3 

He® 

2 ^P Rec 

0.241 

1/4 X 2/3 

He-^ 

1 

1 S Rec 

4.06 

°He++ 

He 

Ly-Oi 

3.0 

2/3 a®^++ 

-J- 

He 

1 1 
2 S - 1 S 

1.5 

g 

2 X 1/3 ajje++ 

He'^ 

Ba-C 

1.0 

B 

®'ne++ 


t Recombinatxon to the 1 S state. 

tt Normalized by n n(X^) . 

e 

- 1 1 

ttt The hh'it of the -IS transition photons which are not 
in the UV have been neglected. 
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the reeombinatxon term, R , is 

B(r) = X [a^(v) + a® aJ 


R.( 

L He 


+ Yy"^ |a^ ^(v) + a® ^ A ^ 


He"^ He'^ 


+ Yy^"^ cP . . ^ 


He 


++ 4*+ 

He He 


.] 


(3.19) 


and the cc terms are approximated by 6 -functions at appropriate 


X 


effective frequencies. For W photons, the A terms are (see Table II ) 


A = 0 
H 




3 A &(v/Vq - IA56) + ^ [^2/3 6(v/vq - 1.522) 


I (0.56) a(v/vQ - 1 . 183 )] 


(3.20) 


A X ^ 6(v/vq ■ 1*5) + 2/3 - 3.0) + §(v/Vq - 


and for OP photons. 

Ah = [2 X I s(v/vq - 0.375) + 2/3 e(v/vQ - 0.75) + §(v/Vq - O.25)] 


A , = V 


-1 


He 


+ 0 


3A 6(v/Vq - O.35I) + ^ [I 6(v/Vq - O. 2 UI) 


1/3 6 (v/v^ - 0 . 285 )] 


(3.21) 


A ^ = 0 
He^ 


In normalized units, the "source function" [Eq. (2.31)] is 



where 


and 



n(r) = n(r^)a^(r) , ^ 


/s 



k-Tt CXg 


( 3 . 22 ) 

(3.23) 

(3-24) 


The value of the dxmensionless parameter , mentioned previously 
in §2.4, determines much of the physics of the model hy specifying the 

A 

relative roles of gas emission and dust extinction. When « 1 , then 
''E/ 

2 (r) IS large and recombination photons dominate dust heating. Alter- 

A 

natively, when » 1 , recombination photons have little effect on 
the dust thermal structure. Once the dust parameters have been specified, 

A 

it IS that determines the total dust extinction optical depth of the 

nebula (see Appendix b). 

The extinction coefficient for UV photons is 


K^(r) = n^(r) t n(r)[(l-x)a^ f^(v) 


(3.25) 


+ Y y®CT ^„«(v) + Yy'^o' . f , (v)3 
He ,0 He ,0 He"^ 

where the radial functions x, y^ , and y^ vary slowly until they reach 


their respective Stromgren radii. In normalized units, Eq. (3-25) is 


K^(t) = f^(v) + [ag(R)/fly(R)] X ijj(v) 

(3.26) 


r 0 

+ 

rid°z-f 

1 e 

„(v) . 1- f ^(v)l 

*- He 

He -J 
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where 


and 


■ V,0 ^ 


He ^0 


'"(^o)‘^o 


(3.27) 


(3.28) 


The parameter e determines the "width" of the hydrogen and helium 
transitxon zones (i.e, regions where the functions x and y go from 
1 to O) by specifying the optical depth of the gas relative to that of 
dust. The functions (i-x)/e j y /e and y /e are all independent of 
e except near their respective Stromgren radii . For normal interstellar 
dust-to-gas mixtures^ e 2 x 10 ^ (and as described before for nebulae 
without dust, e = 1) . 

Equations (3.25) and (3.26) are valid for OP photons if we note 
that f , f , and f are zero except for Ly-a photons. 

H VT ^ TV 

He He 

3 . 7 . Solution of the Ionization Equations 

In dimensionless form, Eq. (3.2) can be written for hydrogen and 
helium as 

x(X^) = x(X^‘^^)[x+Y(y'^ -t-2y’*"^)] A (R)/u(X^) (3.29) 

X^ g 

where co 



+ ^ ^ (3.30) 


4o 



and the summation is over the number, N , of "lines" that can ionize 
the ion X (e.g. for hydrogen, this summation is over all 9 UV lines). 


The coefficients A ^ are 


Ajj = e a. 


B 

0 "H 


B ''0 

A . = e a ^/S- a He 
He He He 


■^w = ® 0! , , /S^ a® . , He'*’ 

He** ° He** 


(3.31) 


Inside of the hydrogen Stromgren radius, 1 - x « 1 and Eq. (3*29) 
is a linear equation for (1 - x)/e . Near the Stromgren radius, where 
X varies rapidly with T , the quadratic form of Eq, (3.29) i-S solved. 
This rapid radial variation of x has been a challenge to all who have 
constructed models of HII regions, Kirkpatrick (I9TT) developed an 
elegant method of handling the variation accurately when the equation of 
radiative transfer is approximated as an initial value problem. However, 
because the exact solution to radiative transfer is a two-point boundary 
value problem, we have developed a less elegant method (in fact, a brute 
force method) in order to improve the accuracy of the ionization structure 
at the hydrogen Stromgren radius ( see §6.4). 

As can be seen from Eqs. (3-29) and (3.30), the function u only 
approaches zero at the Stromgren radius so that x -» 0 as T -* os . It 
IS therefore more convenient and consistent to define the Stromgren radius 


as the point where 



. (3.32) 


* 

tot 




It should be remembered that approaches unity as 0 . It is 

not exactly unity because of the small contribution by J^(T=Oj . 

• 0 

The helium Stromgren radius is defined when y 1 . Inside of 

this point, to avoid solving coupled quadratic equations for y and 

H — h *t" 

y , there are two possible alternatives. First, the y + 2y term 

from n^ can be calculated from the ionization structure in the previous 

0 ■|”I' 

iteration. Second, it can be assumed that y and y are never both 
non-zero at the same radius. It has been found that the first method 
leads to a slightly less rapid rate of convergence of the ionization 
structures. 


3.8. Observables 


A useful quantity in comparing model HI I regions^ with observations 
IS the fraction of ionizing radiation "destroyed" by the gas. The gas 
"destroys" ionizing radiation whenever a recombination takes place that 
does not produce a UV photon. This fraction, ^gas ^ nnder the assumptions 
of §3.1j IS 




gas 


^ ^ o:' 


B 

H 


‘0 


+ Y y'*' X J X 0*^^) a 


B 

He -I 


(3.33) 


where all recombinations to excited states of hydrogen destroy UV protons, 

1 0 / 

whereas only recombinations to the 2 S state of He (where the two 
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photon emxssxon process has a profaabxlity of not producxng a OT 
photon) can destroy UV photons. In. normalxzed units, 


f 

gas 






dR <1 ^(r) 
g 


2 

R 51 
e 


[> 


X + 0.04 Y 


B , B 

a 

He 


]. 


R, 


■0 


(3.3i^) 


and 


= X + Y(y'^ + 2y'*"^) - 


( 3 . 35 ) 


Equatxon (3-3^) can he used to approxxmately detennxne and 

Sq for any spherxcally symmetrxc dusty nebula, xf the radxal varxatxons 

of the gas and dust densxtxes and the spectrum of the central source are 

known. If the dust xs optxcally thxck over most of the spectrum of the 

central source then the total stellar lumxnosxty L xs equal to the 

total observed IR lumxnosxty L , and assumxng a spectrum, S can be 

calculated. If the dust xs optxcally thxn, then xt mxght be possxble 
* 

to observe L and the spectrum dxrectly. Because xt is possxble to 

express approxxmately as an explxcxt functxon of (see Appendxx B), 

then settxng x = y =1 xnsxde the hydrogen Stromgren sphere allows 

Eq. (3*3^) to be expressed as a sxmple function of . In additxon, 

f xs approxxmately proportxonal to the ratio L„^/L „ (x.e. the 
gas ^ X- r IR ^ 

integral xn Eq. (3.33) is proportxonal to the voltime emission 

measure) where the constant of proportionality is the ratio of the 

effective Hp recombination coefficient (AGN p- 66ff) to o:“ . Thus 

Eq. (3.3^) relates an observational quantity to an explxcxt functxon of 
* 

'^1 ' 


* It should be noted that this analysis is not dependent on the determi- 
nation of the gas or dust densities. This xs consistent with our 
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The ratio of the intensity of helium to hydrogen recombination lines 


IS directly observable and is (in the absence of intrinsic redening) 


Rn 


R-, 


■/ 

^0 


L ^ =1 dRa^ (r) R^ny^^/ dRa^(R) r^ji x 


(3.36) 


R 


0 


L . « 1 when the hydrogen and helium Stromgren spheres coincide. For 
+ 

-8 

O stars, L is on the order of 10 

y 4**r 

ti 

Finally, we have defined the average ionization of the ion X as 

R, 


=Jd 


■1 \ 


<x(x ^)> =jdRa‘^(R) r‘^ 71 dRa^ (r) R'^n . (3.3t) 

g 6 IS ® 


R. 


0 


R, 


0 


This quantity is proportional to the volume emission measure of recombi- 

1-1 

nation lines emitted by the ion X 


3 . 9 . Lyman-Alpha Lines 

Radiative transfer of resonantly scattered Lyman-o: lines is 
included in our calculations with a slight modification of Eq, (2.1) and 
subsequently the JAY and RAY equations. Resonance scattering is defined 
to be an absorption (using the doppler broadened absorption cross section 
evaluated at line center) followed immediately by an isotropic emission. 
Equation (2.1) for H Lyman -OJ photons is, for example. 


comments of §3.6 that once the dust parameters are specified, it is the 
combined parameter Sq which determines the physics of the model and 
not the values of the individual parameters which combine to form S 
At this point, however, if the size of the nebula is known, 
can be calculated from Eq, (2.24) once ^2 xs specified. ^ 
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+1 

* ? \^^K,0 V f ’^1^' 

Lj L L 

-1 

+ ( 1 -x) n (r) cr^ (v ) J^ (r) ( 3 . 38 ) 

jj Vl 

where = 0*75 Vq is the frequency of the H Lyman-o: line, 6 ^ is 
along the photon propagation direction, and the last term represents 
the isotropic emission. This equation differs from that of recombination 
lines because the gas source term depends directly on J^(r) . 


45 



4. PHYSICAL PROCESSES IN THE DUST 


4 .1 Basic Concepts 

The parameterization of the dust has been discussed in Section 2.4. 
To reiterate, the dust is parameterized by n^(r), a^(v) and oi and, 
in the case of anisotropic scattering, by and g^. All of these 

quantities are, in fact, averages over a distribution of dust particle 
properties. If the distribution of particle properties is represented 
by a function f(a), a being the particle's radius, then the mean 
radius is 


os CO 



da f(a) 


and, for example, the mean extinction cross section is 




00 


CO 


a^(v) = 






da f(a) 


(4.2) 


0 0 

In a real nebula, f(a) may also be a function of position, time and 
particle composition. Keeping this in mind, we will continue to use 
"mean" particle parameters knowing that these averages can obscure 
some of the physics of the dust. For example, we characterize the 
dust eraissivity by generalized Mie theory so that the emissivity is 


= 3ta (27taA)' 


(4.3) 


where j is in the range 0 ^ j s: 2, and the mean dust radius is in 
the range 0.1|j. ^ a l|j, for normal galactic grains. If we average 
Eq, (4,3) over a particle distribution then an additional factor 
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the order of unity would appear and a would be replaced 
by a . 

For UV and OP photons the frequency dependencies of O'^(v) ^^ud 
03 ^ are given by power laws and for IR photons, both generalized Mie 
theory and a more realistic extinction coefficient from Knacke and 
Thomsen (I 973 ) have been used. 

4.2 Dust Energy Balance and Dust Temperature 

The dust particles absorb UV, OP and IR radiation and radiate 
this energy in the infrared at temperatures from 'v 1000 °K down to 
a theoretical limit of 3°K. Assuming that the heating processes are 
radiative only, the dust energy balance equation is 


GO 

/ 


dvn^(r)ad(v) (1-co ) B ( T^) = 


dyhv n (r)<T (v) (1 -cd )J (r) (4,4) 

a a V \) 


where is the effective temperature of the dust grains at the 

radius r and B (T ) is the Planck function at the temperature T. . 

V d. d 

The left-hand side of Eq. (4.4) represents the infrared emission and, 
in practice, for arbitrary d^(v) the integral extends only over IR 
photon frequencies. The right-hand side calculates the absorption 
of UV, OP and IR radiation. 

To evaluate the dust temperature from Eq . (4.4), it is convenient 
to write the energy absorption term in the form 

Vg 

= VQ^y*dv vfjj(v)(l-cn^)[S^ (t) + R^J®(t)]> (4,5) 

vi 
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where, W (t) is proportioned to the rxght-hand side of Eq. (4.4) and 

jPR 

depending on the limits and slso proportioned to 

the energy absorbed from the UV, OP or IR frequency ranges. In 
actual calculations, (vj^^V 2 ) have the values (vq^^Ovq)j 
( 0.09l2yQ,VQ) and (~j0, 0.0912Vq) for UV, OP and IR photons respec- 
tively. For the lines, Eq. (4.5) integrates to 

w (t) = e\^2v J° (T)(lto )£.(v), (“t.6) 

^ ° j=l ^ ''j '’j “ " 

and IS summed over all I 5 UV and OP lines. F(t) is the sum of the 
W functions : 


f(t) = w^^(t) + Wqp(t) + w^j^(t) + wJt). (4.7) 


The left-hand side of Eq. (4.4) has the form 


/ dx x“^f f x) 

(Tq^T^) > '"0 “ ■"'0''''' == - '-/vo . (•>■&) 

P -1 


f (x) IS the frequency dependence of the infrared emissivity 
IR 

r f_„(v)=(l-co )f,(v)l and T_«15T690 *0°K. Eq. (4.4) now has the form 

'■ IK V) d u 


l(T^) = IF(t)/r" 


(4.9) 


and the dimensionless constant I, which scales the temperature, is 
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A 


(4.10) 





} 


where 


Q. 


gives the efficiency of 


IR emission relative to cr, ^ 

d,0 




(4.11) 


For non-Mie theory emissivities^ is used in Eq. (4.10). Qq is 

the ratio of the extinction to the geometric cross sections and is 

the order of unity for interstellar grains (Spitzer I968, p. 60) . 

As we have seen in Section 3.8, Sq is determined from the luminosity 

of the central source and n,(r )c, ^ can be determined observationally 

a u d,U 

so that Q IS the only free parameter in Eq (4.10) and can be 
varied to scale the dust temperature of the model. In the next 

A 

section, we shall see that I also determines the "strength" of 
the infrared emission. 

For generalized Mie theory, given by 

I(T^) = (V'^0^^^ ^(4tj)£(4tj) (4.12) 


where ^(4 + j) is the gamma function and C(4 + j) is the Riemann 
zeta function. If fjjj(x) has an arbitrary functional form, then 
Eq. (4.8) must be used to tabulate as a function of I. The 

right-hand side of Eq. (4,9) is given by radiative transfer and 
table look-up or Eq. (4.12) are used to calculate T^(t) 
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^l-.3 Infrared Source Term 


The source term of IR photons for the Jay and Ray equatxons is 
given in terms of the IR absorption coefficient by Kirchoff 's law 
(Eqs. 2.2 and 2.3). Using the normalization condition, (Eq. 2.31)^ 
the source function for the Jay and Ray equations is 




(ii-.13) 


The parameter S (with units of 1/Hz), which determines the 
X R 

"strength” of the IR photon emission, is given by 


"iR = 


(4.14) 


So far, S and I are valid for any nebula regardless of the 
XK 

existence of ga^. For nebulae with gas more conveniently 

expressed in terms of the gas density and dust-to-gas ratio as 




Ho 




Rj3Gfl0^c£3/n(r^)] [2«/0.9l2]' 


^0^ 


(^.15) 


0 


M is the mass of a hydrogen atom and the "dust-to-gas" ratio, R „ , 
■d DG 

IS defined as 


R^_ = 
DG 


M M 
d g 

0.01 


Igr/cm' 


O.lja 


J 


J-1 


(4.16) 
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where M is the dust-to-gas mass ratio and n is the density of 

d g 

o 

a grain (p « 1 gr/cm"^) . For non-Mie theory, Eqs. (4.15) and {h ±6) 

A 

are evaluated for j = 0. With S and I m the form of Eq, (4.15), 

IK 

it is the parameter which can be freely adjusted. Combining 

DG 

Eqs. (3.28), (4.11) and (4.l6), it is easy to see that 


e = (2 X 10'^(2rt/0.9l2)^Rjj^/Q^ . (i^.l?) 

A 

From Eqs. (4.9) and (4.13) it can be seen that S (or equi- 

IR 

A 

valently l) determines not only the "strength" of the IR emission 

but also the spectrum. By noting that neither Q nor R„_ depend 

j JUG 

on Sq, it IS easy to see that the variation of these parameters will 

change the qualitative shape of the emitted infrared spectrum while 

keeping the total IR luminosity, constant. We will return to 

IR 

this point in the next section. 

The observed flux of infrared photons is 


A p A T> ^ ^ 

= FE (H^ t ) 




where 


F = hvQSQ/4«D^ . 


(^.19) 


D IS the distance to the nebula. S„„(v) is usually expressed in 

IK 

-23 2 

Janskys (JY) where one Jansky is 10 ergs/cm -sec-Hz. is 

obtained from the total luminosity of the source which, in turn, is 

2 

obtained from a measured luminosity times 4rtD . Thus the quantity 
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A 

Sq/D xs usually independent of distance. I is also independent of 
distance if q is calculated from an observed angular size. 

Generalized Mie Theory 

For many HII regions that are heated by 0 or B stars, the 

dominant source of dust heating is from UV and OP photons. Also, 

in many nebulae the infrared spectrum is much broader in frequency 

than a single black body indicating that the dust must not be optically 

thick over most of the spectrum, if t » 1 for IR frequencies, then 

> V 

we would only see into the nebula to a depth where t 1 a^ud the 

V 

spectrum would necessarily be much narrower. In these cases it may 
be physically reasonable to assume that the dust is optically thin 
to all IR photons. This assumption is very useful for the purposes 
of understanding this model because it allows the investigation of 
the effects of radiative transfer in the UV and OP frequency ranges 
on the dust temperature structure while preserving the quantitative 
features of the IR spectrum 


When the dust is optically thin to IR photons, the observed 
spectrum is 


■/' 




(4.20) 


0 


which in terms of S and F is 

XR 




^IR^ "O' 


T- 

/; 


djR 


ThvTkO 


d' -1 


(4.21) 


0 


T (r) IS calculated from Eqs. (4.7), (4.9) and (4.12) with W = 0. 
u XR 
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The total xnfrared lurainosxty of the source xs 


/■ 

0 


^IR "" I ®rr>(v) . 


IR^ 


(4.22) 


Suhstxtutxng Eq. (4,21) xnto (4.22) and changxng the order of xntegra- 


txon, L IS equal to 
IR 


^IR ” ^ ' dT ^(t), 


T- 

f 


(4.23) 


independent of S (or equxvalently, of l). It can be explxcxtly 

IK 

A A 

seen that varxatxon of S and I through the parameters Q or 

iK j 


R_„ does not change L • 

iJCj IR 


In addxtxoHj whxle S affects the 

IR 


qualxtatxve shape of the IR spectrum, xts value cannot be obtaxned 

observatxonally. Thus Q or R can only be obtaxned by "fxttxng" 

J DG 

an observatxonal IR spectrum wxth the results of the model. 

Eq. (4.21) assumes that the nebulae just fxlls the beam of the 
observxng telescope. In general, however, the flux xs measured as a 
functxon of aperture sxze where the aperture of the telescope xs less 
than the sxze of the object. In the optxcally thxn case, the IR flux 
from wxthxn a cxrcle of radxus b, centered on the nebula, xs 


IR' 



;r^-t^)"2 j (r) b > r 
V 0 


(4.24) 


b :£ r^ 
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Changing the order of integration, Eq. (4.24) reduces to 


4jtD S (b,v) = 4rt 

IK 


X 

[/ 


dr r j (r) - / dr r 


Pk ^ 

■v/r - b 


J (r) 
V 


^0 


maxCb^r^) 


where max(b,rQ) picks the larger of (b^r^). When b = r^, 
Eq. ( 4 . 25 ) reduces to Eq. (4.20). 
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5. APPROXIMATE SOLUTIONS OF RADIATIVE TRANSFER 


5 . 1 . Introductxon 

Approximate solutions of radiative transfer have been used exclusively 
in previous models of the ionization structure of H II regions. Hummer 
and Seaton ( 1963 ) were the first workers to calculate accurate models 
using the on-the-spot approximation (osa) for planetary nebulae (without 
dust) . Tarter (I967) developed a more elegant solution using the "outward 
only" approximation. For models of H II regions with dust, Balick (1975) 
has used the OSA and Sarazin (19TT) used the "outward only" approxi- 
mation. All of these solutions are unable to handle dust scattering 
explicitly. In this chapter we will discuss approximate solutions of 
radiative transfer in some detail including the OSA, the diffusion 
approximation (valid for nebulae with little or no gas), and a "modified 
Eddington approximation". In Chapter 7^ we will make detailed compari- 
sons of the results of these approximate solutions with the QDM solution 
in order to assess the accuracy of the approximations. 

Petrosian and Dana (1975^ Paper l) and Dana and Petrosian (197^^ 

Paper II) have developed semi -analytic approximate solutions of radiative 
transfer y/hich are accurate to within 50 ^^ which are convenient for 

I 

the initial modeling of dusty nebulae. The important results of these 
papers will be presented here. 

5.2. On-the-Spot Approximation 

The existence of a Stroragren sphere requires that all ionizing 
photons be absorbed in the nebulae, i.e. the number of diffuse ionizing 
photons emitted per second must equal the number of diffuse photons 
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absorbed per second. Thus, for UV photons. 




/• 


S' / \ D/ \ 
dVK^(r) J^(r) 


(5.1) 


where the integrals are over the entire volume of the nebula. The 
essence of the on-the-spot approximation is to assume that, in the 
absence of dust scattering, Eq. (5-1) is valid locally so that 

j^(r) = ic“(r) J®(r) . (5.2) 


We shall distinguish between two different forms of Eq. (5.2). First, 
the generalized on-the-spot approximation (gOSA) assumes that the diffuse 
intensity is absorbed by both gas and dust. Second, the on-the-spot 
approximation (osa) assumes that the diffuse intensity is absorbed by the 
gas only. 

For recombination photons to the ground state, the effective frequency 
is very close to the ionization limit of H^, He*^, and He^ where the 
respective photoionization cross sections are the largest. The gas 
opacity IS therefore large for these photons and they should be absorbed 
close to where they were emitted. This is also true for Lyman-a photons 
which are resonantly scattered and hence cannot travel far before they 
are absorbed by dust. Equation (5*2) is not valid, however, if dust 
scatters UV photons or if the nebula is density bound (i.e. there is 
insufficient gas or dust in the nebula to absorb all UV photons) . 

5 . 2 . 1 . Generalized On-the-Spot Approximation 

From Eq. ( 2 .I 5 ) with = 0 , it can be seen that Eq. (5*2) 

D/ \ \ 

implies that S (T) = 0 , and that the normalized mean intensity, J. (.'T; , 

S) V 
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xs given by the ratio of Eq. (3-22) to Eq. (3.26). Wxth given by 

GOSA and given by Eqs. (2.45) and (2.49)^ the ionization structure 

and dust temperature functions W^(t) and Wj^('^) ^^.n be calculated. 

Other solutions of radiative transfer must be used, however, for OP and 
IR photons. A GOSA solution has been included in our numerical model 
thq^t has the same gas and dust properties as the QDM solution. 

5.2.2. Semi -Analytic Approximate Solution 

We have developed semi -analytic OSA models of dusty H II 
regions containing hydrogen and helium. These models use frequency and 
angle averaged cross sections and intensities. Approximate analytic expres- 
sions are derived for evaluation of the variation of the fraction of lonizii 
radiation absorbed by gas and the ratio of the volume emission measures 
of He II to H II regions with the spectrum of the central source, helium 
abundance, and absorption properties of dust. For nebulae containing 
hydrogen and helium, we have derived simple coupled differential equations 
for the net flux of ionizing photons and for the fraction of ionizing 
photons capable of ionizing helium. These initial value problems are easily 

A 

solved numerically to give S (t) [Eq. (2.16)] and the ionization structures 
of hydrogen and helium. In the absence of helium or in the absence of 
dust, these equations have analytic solutions. Appendix B gives the 
derivations and important results of the OSA models. 

5 . 3 . Diffusion Approximation 

The diffusion approximation has often been used in models where 

scattering dominates radiative transfer. The approximation assumes a 

single isotropic point source surrounded by an isotropic medium in which 
s ^ 

K: » K and the scattering is isotropic. The medium can either be 

V 
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infxnite or be surrounded by a "free surface" x.e. a surface on which no 
photons enter from outside. The solution of the diffusion approximation 
IS valid a few scattering mean free paths away from the central source 
and away from the free surface. 

The inner boundary condition is 

S^(t=0) = S* ( 5 . 3 ) 

and for an infinite medium, the outer boundary condition is 

Limit J (r) = 0 . (infinite Medium) (5*^) 

r eo 


At a free surface, it is assumed that vanishes at an extrapolation 

distance d (Davison 1957^ §8*5) outside of the surface where d is 
obtained from the exact solution, of Milne's prgblem (infinite half -space) . 
The free surface boundary condition is 


J^(r^+d) = 0 , d « 0 . 7104 /K^ , (Free Surface) (5*5) 

where 1/*^^ the extinction mean free path. The key assumption in 
this approxim.ation is that d/r^ « 1 . 

We will assume a uniform medium. However, this assumption is not 
ma_thematically necessary and nonuniform nebulae will have a different 
numerical factor in the "diffusion constant" which comes from the value of 

r 



(5.6) 


i ^ 



) 


integrated oyer the volume of the nebula. 
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The diffusion equation for J (r) is 

V 

J^{r) , =Jy^ . (5.T) 

S S' 

K and K are functions of frequency hut are assumed uniform throughout 
the medium. In spherical coordinates, the solution to Eq. (5*7) 

Jy(^) = i (Ae ® + Be^ . (5.8) 

The ratio A/B can he obtained from the outer boundary conditions on 
J^(r) [Eqs. (5-^) or (5-5)] • I'll® n®* flux through a spherical shell is 
given by Fick's law 


S^(r) = - (Ij-jrr^ 73^^) r*V[ivn J (r)] 


( 5 . 9 ) 


Equation (5«9) is derived assuming that the medium is infinite (i.e. it 
IS not valid within a few mean— f ree-paths of surfaces) and assuming that 
J^(r) is a slowly varying function of position. A nice derivation of 
Fick's law IS given in Lamarch ( 1966 ^ §5*^)* Upon substituting Eq. ( 5 . 8 ) 
into (5-9) ^ud using the inner boundary condition [Eq. (5*3)]^ the flux is 


S^(r) = S* - e(l-Tj,)e’'®] /(1-e) , 


( 5 . 10 ) 


where 


e = 


jexp ^-2 \/3ui3(l-^J5) (t^ + 0.71o4)J 


0 


(Free Surface) 
(infinite Medium) 

( 5 - 11 ) 


If T is frequency independent, then S is replaced by S 


0 ' 
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5.^. Modified Eddington Approximatxon (mEA) 

Hummer and Rybicki (I 97 I) introduced a modified Eddington factor 


which they used as an initial estimate of f (r) [Eq, (2.10)]. 
approxiraationj we solve the JAY equation with 

-T 


In this 


f (r) = 

V 


1 -(2/3)e 

1/3 


Continuum frequencies 
Lines 


(5.12) 


for all iterations. This approximation has the correct form for the 

continuum f (r) at both the inner boundary where is isotropic 

V ^ 

and at the outer boundary where is radially outward (f 1) . 

V V 

Because the lines are emitted isotropically^ their intensity remains 
fairly isotropic throughout the nebula. 
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6. COMMENTS ON PROGRAM SCATER 


6 . 1 . Introduction 

Program SCATER is the numerical model that has heen constructed 
to solve the problem of radiative transfer in dusty nebulae. The program 
currently resides in the NASA-Ames Research Center’s CDC76OO computer 
and IS available for public use. SCATER is currently about 6OOO cards 
long and requires 25; 000 words of small core memory and ~ 75, 000 words 
of large core memory. ^Vhile SCATER was being written, an effort was 
made to make the code as efficient as possible (within the author's 
limited knowledge of code optimization) and most Eddington approximation 
models can be run with less than 30 seconds of CPU time. However, because 
the program was written for a general dusty nebula with or without gas, 
and because it was written over a period of about a year, there are many 
parts of the program which are not necessary and/or efficient when the 
program is applied to a specific nebula. Also, the program is not complete 
as it does not calculate a self-consistent electron temperature or the 
forbidden lime emission from the gas. It is estimated that this additional 
physics would require 2 to 3 months of additional programming and testing. 
Inasmuch as this thesis is a users' manual and guide to SCATER, we will 
deal briefly in this chapter with the iteration procedure, timing, and 
accuracy of the numerical model, 

6 . 2 . Iteration Scheme 

Figure 3 gives a brief flow chart of the iteration scheme of program 
SCATER. Most of the "bookkeeping" details are necessary because the 
program uses the same subroutines for UV, OP, and IR photons.,. This 
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Fig. 3. Flow Chart of the Iteration Scheme of Program Scater 
Circled numbers refer to comments in the text 
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iteration procedure is necessary because the source functions 2 ^(t) 

^ s 

[Eq. ( 3 . 22 )] and 2^ (R^b) [Eq. (2.43)] depend on the total intensity 
and must be updated each iteration. In addition, if there is gas, the 
extinction coefficient is also dependent on the total intensity and the 
iteration procedure is much less rapidly converging. 

The following numbered comments refer to the numbered sections of 
Fig. 3 . 

1. On the first iteration, the Eddington approximation or the 

modified Eddington approximation [Eq. ( 5 . 12)] is used to evaluate f^(T) 
[Eq. ( 2 . 10 )], [Eq. ( 2 . 11 )] and f® [Eq. (2.52)]. If there is gas, 

then the initial hydrogen ionization is given by the OSA analytic model 
(see Appendix b) using and the absorption optical depth. If 

CD ^ 0.5 for UV photons, then using cd^ = 0.5 for the absorption optical 
depth in the OSA model produces a better guess for the value of at 

the Stromgren radius. 

2. Once the initial ionization structure has been calculated, the 
stellar flux as a function of T and v can be calculated. If there 
IS no gas or if OP or IR photons are being considered then the stellar 
flux does not depend on the ionization structure (extinction is by dust 
only) and is calculated only on the first iteration. 

3 . The UV and OP recombination lines source term is calculated 


from the ionization structure, and the JAY equation can be solved for 

jJ(T) , T ^ . 

i .t- 

4. For nebulae with gas and for UV photons, the Stromgren radius 


condition [Eq. (3.32)] on the flux is checked at the old value of . 
If ^tot^**"!^ large, an extend flag is set (meaning that the 
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Stromgren sphere will be moved outward) and ) is calculated for 

T > assuming radially outward flow and extinction by dust only. 

5 . The new value of used to update the scattering source 

function . Using ^(t) calculated in step 3, the RAY equation 

IS solved within the old Stromgren sphere. 

6 . If the extend flag is set, U^(r,b) [Eq. (2.38)] is extended 

beyond the Stromgren sphere assuming dust extinction and no sources for 
R > IS the dust extinction optical depth along a ray starting 

at the old Stromgren radius. Eor B > R^ , it is assumed that U^(R,B) = 0 

7. Because IS calculated from the derivative of U^(R, b) 

[Eqs. (2.40) and (2.47)], 'the program avoids its calculation unless it is 

required in the scattering source function (i.e. 7^ O) or it is the 

^ 'H' 

last iteration and the diffuse flux is desired for comparison with S^(t) . 

8a. At this point, frequency integrals are calculated assuming that 
the integrands vary as power-laws between two adjacent frequency grid 
points. The total integrated mean intensity, J(t) p is calculated in 
each frequency bin and is equal to 


m 


j(t) = f dv [s*(t)/R^ + J^(T)1 (T) , (6.1) 

V-, 


where ra is the number of lines in the bin and is equal to zero for IR 

photons or for nebulae without gas. 

8b. For nebulae with gas and for UV photons, ^tot^”*^^ calculated. 

A new value of is defined as the first value of T for which 

S < 10 . The ionization functions for hydrogen and helium are 

tot 

calculated for T ^ [Eq. (3.30)] 
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8c. The- dust energy fun<?t3Lons, W (t) [Eqs. (^. 5 ) (^*6)]. are 

j R 

evaluated on the last iteratxon for UV and OP photons. If the dust xs 
not optically thin to infrared photons, then dust tempera- 

ture are calculated on each iteration of the IR intensity. 

9 . J(T) IS compared with its values on the previous iteration. For 
UV photons and nebulae with gas, 1-x [Eq. (3.8)] is also compared with 
its previous values. If both of these functions have converged to within 
a specified error limit throughout the nebula, the program does one more 
iteration. 

6,3. Timing and Array Sizes 

The timing of program SCATER is critically important if it is going 
to be used to fit the spectra of observed objects. The fitting process 
requires many runs as several parameters are usually varied. In general, 
the ODM solutions take too much times to be of value in the fitting 
process. However^ either the Eddington approximation or the modified 
Eddington approximation run in a reasonable amount of time and are 
accurate (see Chapter 7) • 

Because of the large variation in optical depth scales associated 

2 

with various frequencies and because of the 1/r dependence of the 

diffuse mean intensity at the center of the nebulae, it is not practical 

to use an optical depth distance scale. We therefore used a uniform 

distribution of radial grid points, with extra radial points added at the 

2 

center of the nebulae to handle the 1/r dependence and for models with 
gas, extra radial grid points are also added at the Strbmgren radius. We 
used loo radial grid points and llO impact parameters . Each frequency 

* Dr. Alan Tokunaga of NASA-Ames has doubled the size of the radial 
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band has up to kO grid points. In the continuum, the frequency grid 
points are determined by spectrum of the central sour_ce and by the edges 
of the gas. 


In order to make useful comparisons, a standard model without dust 
scattering has been run to calculate the ionization structure of hydrogen. 
Without scattering or helium, only one frequency grid point is calculated 

p 

for UV diffuse photons (i.e. the HI S recombination line). Each model 
calculated only W photons and executed 30 iterations. The times of 
execution (not including compiling times) were 8.638 seconds for QDM, 

3.264 seconds for the Eddington approximation, and 3.256 seconds for the 
GOSA. The QDM time can be reduced to 6. 815 seconds if the Eddington 
approximation is used for the first I 5 iterations while the program 
"finds" the Stromgren radius. The difference in time between the QDM 
and the Eddington approximation is approximately proportional to 
1/2 X frequencies) , i.e. the time required to solve for on 

the (R,b) grid for each frequency. Additional time is also required to 
calculated J^(t) and f^(T) using Eqs. (2.46) and (2.48). Thus, if 
scattering (tD^y = l.O) is included, the QDM time increases to I 87.7 
seconds while the Eddington approximation time is only 7.964 seconds. 

6.4. Stromgren Sphere 

The rapid radial variation of the ionization structure at the 
Stromgren radius causes serious problems of the accuracy and rate of 

^ (Contd.) 

arrays in SCATER as part of his effort to use the program to fit the 
infrared spectrum of Sl‘40. 
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convergence of the numericaX model. Interior to the Stromgren radius, 
the ioijization structure is rapidly converging and accurate. To deal 
with these problems, extra radial grid points are added where the hydrogen 
ionization varies most rapidly. However, because the exact position of 
the Stromgren radius is not known a prion, these points are added several 
iterations into the TJV calculation when the value of is no longer 

changing . 

The new radial grid points are added symmetrically around the grid 
point where the variation of 1-x from the previous iteration is maximum. 
Depending on the number of unused grid points in the radial arrays, 2, 8, 
l8, or 32 points are added between the first through fourth grids, respec- 
tively away from the point of maximum variation. The functions 
f^Cx) and H^(t) are estimated at the new grid points using a 
cubic splines fit through the old points. The ionization functions are 
estimated assuming a power-law dependence between the old grid points. 

A better "fitting" function could be used that might speed up the conver- 
gences of 1-x . Once this procedure is completed, the B grid is 
adjusted so the B grid is equal to the R grid and a new iteration is 
started. 

6.5* Accuracy 

During the construction of program SCATER, each subroutine was 
carefully tested and comparisons were made with analytic results whenever 
possible. In addition, there are two models that can be run to test the 
differential equation solvers and dust physics and to test the gas physics. 

For pure dust nebulae with = 1 , the total net flux through a 
spherical shell must remain constant throughout the nebula. To within 
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1 to 2^^ over most of the nebula the model kept the fl\ix constant where 

the difference is attributed to the linear approximation of b) 

between B grid values. This effect can be corrected (e.g. Mathis 1972) 

but the error is well within acceptable astrophysical error limits. 

Tor pure hydrogen nebulae, direct comparison can be made to Hummer 

and Seaton's (I 963 ) on-the-spot calculations. Using our built-in GOSA 

model, we find close agreement between the hydrogen ionization, 1-x , and 

the optical depth of hydrogen at the Lyman limit, except at the Stromgren 

radius. Because Hummer and Seaton take into account the frequency dependence 
2 

of H 1 S recombination spectrum and because of differing radial grids, 

close agreement is not expected at the Stromgren radius where diffuse 

photons are the dominant source of ionization. In addition, our model 

calculates f [Eq. (3-33)] to be ~ 2'^ below the exact value of unity 
gas 

in this case. This error is attributed primarily to the rapid variation 
of Ug(r) [Eq. (3-l8)] at the Stromgren radius. 
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7 VARIATION OF PARAMETERS AND 
COMPARISONS mTH OTHER CALCULATIONS 

7.1 Introduction 

Although we have calculated model dusty nebulae usxng program 
SCATER for a mde range of the parameters, we will limit this dis- 
cussion to those questions which have not been addressed in previous 
calculations. 1 ) IVhat are the effects of dust scattering on the 
internal structure and observable parameters of dusty nebulae‘s 
2 ) How accurate are approximate solutions of radiative transfer”? 

The internal structure of gaseous nebulae without dust is well under- 
stood theoretically from the work of Hummer and Seaton (1963^196^)^ 
Hjellming (I966) and Rubin (I968) . More recently, Balick ( 1975 ) 
Sarazin ( 197 T) iis-ve analyzed in great detail the effects of dust 
absorption on the ionization and thermal structure of dusty HII 
regions. Neither of the calculations, however, include scattering. 
Approximate models of dusty HII regions that analyze the effects of 
UV absorbing dust have been calculated by Mathis (I97O, 197!)^ 
Petrosian, Silk and Field (1972)^ Krugel ( 1975 ) j Natta and Panagia 
(1976) and by Papers I and II. 

In this discussion we will consider two types of models. Type I 
models have a constant Stromgren radius while the dust parameters are 

A 

varied. In these models, S^ (in practice, Sq) is "tuned" to keep the 
Stromgren radius fixed. Type II models have a constant value of Sq 

A. 

(or Sq) while the dust parameters are varied. 

We define our standard model to consist of a uniform nebula at a 
distance of 1 kpc with a constant dust-to-gas ratio. The stellar 
spectrum is a i)- 0 , 000 °K (log g = ^l.O) model atmosphere from Auer and 
Mihalas (I972) . Once the dust parameters have been specified, Sq 



determxnes the properties of the gas (see Section 3.6). We have 

A 

chosen Sq = 4.13 which corresponds to (hut is not uniquely defined 

by) a single zero-age main sequence central source with radius given by 

3 - 3 -4 

Panagia (1973) , hydrogen number density of lO'^cm , and e x lO” • IVhen 

A 

Y = 0.1, and ^ photons, this Sq gives 

T*! 1.5. 

Table III summarizes the physical parameters that must be specified 
to completely determine the properties of our model. If the model 
nebula contains no ionized gas, then the parameters S^, Y, e and 
G (r) do not need to be specified. Conversely, if the model nebula 

A 

contains no dust, then the parameters f (v)^ (d } g-i ^ and s do 
not need to be specified. However, e = 1 and G^(r) = C (r) in 
this case. The values of the parameters for our standard model 
(model A in the next section) are also indicated. The value of 
IS entirely arbitrary and was picked for numerical convenience. 

7-2 Variation of Albedo and Extinction Coefficient 

We will consider the effects of the variation of the albedo and of 
radiative transfer on the mean intensity and Eddington factor using 
type II models (constant central source luminosity). The same 
quantitative effects that we will discuss are also seen in type I 
models (constant Tj^) . 

Also in this section we will discuss the effects of the variation 
of the albedo and the frequency dependence of the extinction coeffi- 
cient on the ionization structures of the gas for type I models. The 
observational question we are trying to answer is If I have line 
intensities from a resolved HII region, how do dust parameters affect 
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CQ> cttjcnjS Hi M> tdQP5t<m m> 


TABLE III 


PHYSICAL PARAMETERS 


Parameter 


D 





L 

2 


+++ 

IR 


Defined 

Model A 

Comments 

3.24'^ 

4,13 

Determines ^.nd "relative strength" 

of dust and gas physips 

3.28 

_4 

2 X 10 

Determines "width" of Stromgren radius 

§3.1 

0.1 

n(He)/n(H), He abundance relative to H 

2.30 

1.0 

Radial variation of dust and 

3.23 

1.0 

gas densities 

2.24 

0.1 

Radius of central core in optical depth 
units 

§2.4 

* -H- 

T = 

AM 

40,000 

Normalized spectrum of the central 
source 

2.25 

1.0 

Frequency dependence of the dust 
extinction coefficient 

2.26 

0.0 

Albedo of the dust 

2.6 

0.0 

Average cosine of the scattering angle 

2.7 

0.0 

Average of the square of the cosine 
of the scattering angle 

4.15 

0,356 

Determines the dust temperature and 
the shape of the IR spectrum 


I 


I 


+ Equatxon number 


++ 


Auer and Mihalas (1972) Spectrum, 
photons cm sec , r^ = 5*^^ x 10 cm 


r^ = 1.667 X 10 

(panagia 1973) 


24 


UV 


+++ For j = 0 
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the physical parameters that I deduce from ray observations'? 

Finally, we will discuss the effects of the variation of the 
albedo on the observed IR spectrum again using type II models for 
clarity of interpretation. However, it is usually the infrared 
luminosity of a source that is known> observationally rather than 
the central source luminosity. Hence, these results are only a 
qualitative guide to the effects of albedo. 

7.2.1 Mean Intensity and Eddington Factor 

On Figure 4 , we have plotted the variation of log [ J^(t)/J^( 0 ) ] 

versus t / tj _ ~ '^0 ^ "tliese photons can photoionize hydrogen) , 

y = Vq (these photons cannot photoionize hydrogen), v = 0.25vQ(Ba-C) 

and \) = 1 . 06 vq (H l^S recombination "line") . The difference between the 

shapes of the continuum curves (v = Vq) and the line curves 

(y = 1 ,o6vq) is due to the difference in the "source functions" of the 

two types of photons. The source of continuum photons is, for isotropic 

scattering, 1C®J ( t) which is a strongly decreasing function of t whereas 
V V 

the source of lines is, for hydrogen, xnnJDJjj which is almost constant 
inside the Stromgren sphere for uniform densities. 

For UV photons (y = vq and y = 1 . 06 yQ), it can be seen from the 
figure that as albedo increases the relative intensity is lower near 
the center of the nebula and higher as t approaches Thus the 

effect of albedo is to "transport" photons from the center of the 
nebula to the outer regipj^.. This is primarily a result of the fact 
that we are keeping the extinction coefficient constant while the albedo 
is increased and we are, therefore, also increasing the dust absorption 
mean free path. For type I models if we increase the albedo while 
keeping all other parameters fixed, J^(t) will increase at the same 
radius because the absorption coefficient is decreasing. This effect is 
most apparent in the value of (i.e.'the size of the Stromgren sphere). 
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T/T, 

Fig. 4. Variation of the Mean Intensity Versus Optical Depth. 

Variation of Log [J (t)/J (0)J co - 0 (solid lines), Oi = 1 
(dashed lines) and ^for (H loHizjng photons), y = '^y" 

(non-ionizing photons), y = l,o6yQ (H 1 o recombination "line") 
and y = 0.25 vq (Balmer continuum ^line”) . Circles indicate 
the yadii where Ttt = 1. 

H 
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For these type II models with Y = ) and 8^=4. 13^ co_^= 0 and 

T-, = 2.388 for oj =1. Hence for the same central source luminosity^ 

-L V 

variation of the UV albedo from 0 to 1 makes the volume of the 
Stromgren sphere 2.75 times larger. 

When the albedo is zero, only stellar photons contribute to the 
mean intensity at v = Vq • Thus the difference between the two 
curves is due to gas absorption of \j — Vq photons . On the = Vq 
curves we have indicated the points in the nebula where the optical 
depth of neutral hydrogen at V = Vq, equal to unity. Tjj 

IS defined as 



r 

y (1 - x) n(r') CT^dr' 


(7.1) 


Interior to that point, the ^ curves converge because the gas 

3 

IS optically thin. At the Stromgren radius, t, rv 10“^ in all of our 

H 

models . 

Because the gas does not absorb V =* Vq photons, the large 

difference between the curves for a, = 0 and 1 is due 

entirely to the properties of the dust. When ca = 0, (r) <= e”*'' 

V Vq 


where t is the dust extinction (and absorption) optical depth. How- 

o 2 D 

ever, when = 1, r h (r) = r(H +H ) = constant. Thus we 
V Vq Vq Vq 

see that the effects of albedo on continuum mean intensities are much 

greater for photons that are not absorbed by gas. In addition, the 

"source function" for diffuse continuum photons is proportional to 

(JdJ (t) a-nd hence the mean intensity is strongly dependent on the 
V 

albedo. This is not true, however, for lines as can be seen by 
comparing the y = 0.25vq curves for oi = 0 and o) = 1, In this 
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case, the varxatxon of the mean xntensxty xs almost xndependent of 
albedo because the "source functxon" xs roughly constant. 

Although we have npt plotted the varxatxon of the mean xntensxty 
for "hard" UV photons (v > ^Vq)^ should mentxon here that the 
xonxzatxon at the Stromgren radxus xs sustaxned by these penetratxng 
photons (whose gas opacxty xs small over most of the nebula) and by 
recombxnatxon photons . 

The Eddxngton factor, calculated self-consxstently by the QDM, 

depends on the boundary condxtxons and on the source of the photons. 

Our boundary condxtxons assume that there xs no source of photons 

from outsxde of the nebula. Thus dxffuse contxnuura photons must all 

be travelxng "outward only" near the outer boundary and f (<r) 1. 

V 

However, for Ixnes, whxch are emxtted xsotropxcally and whose 

emxssxon does not depend dxrectly on (t)j f (’t) stays close to 

the xsotropxc value of ^ throughout the nebula. On Pxgure 5^ we 

present the Eddxngton factor from the standard model for the dxffuse 

contxnuum at the Lyman ixmxt and for hydrogen n = 1 and Balmer 

contxnuum recombxnatxon Ixnes. It can be seen that even for xnter- 

medxate values of and <r (t ^ 1.5 for the standard model), the 

V 

dxffuse contxnuum becomes sharply forward peaked at the outer boundary, 
and the Eddxngton approxxmatxon should dxffer from the QEM solutxon xn 
thxs region. The modxfxed Eddxngton approxxmatxon (see sectxon 5*^) 
does a better job of approxxraating the varxatxon of the Eddxngton 
factor for both the contxnuum and the Ixnes. 

Physxcally, the reason that UV contxnuum photons become so 
forward peaked at xs that those photons whxch are not travelxng 
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radxally outward have a longer path length in the gas and are hence 
preferentially absorbed by the gas. It should be remembered, however, 
from Figure 4 that the intensity of these photons is so low at 
compared to UV lines that they do not contribute to either gas ioniza- 
tion or dust heating. Because lines have a constant source function, 
they can maintain isotropy throughout most of the nebula. The OP line, 
which IS not absorbed by gas, does "see" the edge of the nebula and 
becomes somewhat forward peaked right at the Stromgren radius. On the 
other hand, because of the large gas optical depth, the v = 
line does not "see" the boundary and remains nearly isotropic through- 
out. The Lyman-o: line, for which we have not plotted the Eddington 
factor, maintains f (t) = throughout the nebula. 

Finally we note that the variation of the Eddington factor is 
almost indepedent of 
7.2.2 Ionization Structures 

The overall ionization structures of hydrogen, helium and the 
trace elements are presented in Table IV. These type I models, 
calculated using the Eddington approximation, have the same central 
source spectrum and Stromgren radius while the albedo and frequency 
dependence of the dust extinction coefficient (Eq 2 25) are varied. 

We have used a power- law dependence for f^j(v) 

fjj(v) = (v/Vq)^ (7-2) 

where 5 is constant for both UV and OP photons Models A, B and C 
have (0^ s 0 and 6 == 0 , 1 and 2 respectively. Model D is a nebula 


77 




TABLE IV 


IONIZATION STRUCTURE OF MODEL NEBULAE 


Models 


Parameter 


B 

c 

D 

E 

F 

/s 

®0 

4.128 

9.0 

13 0 

1.69x10^ 

2.5 

1.3 


0 

1 

2 

0 

0 

0 

C0 

V 

0 

0 

0 

0 

0.5 

1.0 

f 

gas 

0.298 

0.158 

0.102 

0.927 

0.495 

0.961 


1.015 

0.889 

0.383 

1.019 

1,015 

1.014 

<x> 

0.973 

0.995 

0.991 

0.977 

0.974 

0.976 


0.988 

0.884 

0.380 

0 995 

0.989 

0.989 

(o"> 

0.035 

0.254 

0.793 

0 036 

0.030 

0.027 

OJ 

0 

0.952 

0.740 

0.196 

0.953 

0.957 

0.962 

(c"") 

0.o43 

0.157 

0.506 

0.021 

0.039 

0.036 


0.426 

0.642 

o.ii66 

0.451 

0.398 

0.375 

(c+3> 

0.525 

0.199 

0 024 

0.505 

0.558 

0.586 

<N^ > 

0.o43 

0.224 

0.691 

0.032 

0.038 

0.035 

CM 

0 818 

0.719 

0.286 

0.845 

0.818 

0.819 

<n*3) 

0.135 

0.054 

0.0091 

0 107 

o.l 4 o 

0.142 

(s'") 

0.0097 

0.036 

0.155 

0.0021 

0.0090 

0.0083 


0.068 

0.314 

0,680 

0.069 

0.059 

0.051 

<s*3) 

0 733 

0.577 

0.150 

0.767 

0 . 73 ^^ 

0.735 

<s•^^ 

0.188 

0.071 

0.011 

0.156 

0.198 

0 205 

(Ne"^ > 

0.133 

0.479 

0.899 

0.136 

0.117 

0.104 

(Ne 

0.864 

0.516 

0.083 

0.864 

0.881 

0.893 
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TABLE IV (Cent.) 


Models 


Parameter 

A 

B 

c 

D 

E 

F 

(Ar^> 

0.0056 

0.117 

0.612 

0.0017 

0 .0051 

0.0046 

<Ar+=> 

0 294 

0.543 

0.315 

0.306 

0.266 

0.243 

<a/^> 

0.698 

0.337 

0.052 

0.692 

0.727 

0.750 
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wxthout dust. Models E and F have 6=0 and oi =0.5 and 1.0 

V 

(xndependent of frequency) respectxvely . The functxons f apd L 


gas 


+ 


He 


are calculated with Eqs. (3.3^) through (3*36) and xn all of these models. 


Xj ^ 

RJ 10 . The average xonizatxons are calculated using Eq. (3.3T) 

and we have not listed ions for which (X^^) < lO"^. 

In models A and D through F, L ^^1.0, in agreement with the 

He 

results of Hjellming (I 966 ) but not with Hummer and Seaton (I 964 ). 

Sarazin (l9TT)^ who finds L ^ < 1.0, incorrectly concludes that this 

He 

IS an effect of the on-the-spot approximation. Our parametric studies 

have shown that the details of the ionization structure at the 

Stromgren radius depend most strongly on the spectrum of J ( t « Tt ) • 

V 1 

However, as Sarazin states, the question is academic because static 

photoionization models are not valid at the Stromgren radius where an 

ionization front and associate shock waves exist. We note also that 

L ^ does not depend on cd • Because to ;£ this conclusion is 
He V V 

independent of the frequency dependence of the albedo. As 8 
increases dust competes more effectively for He ionizing photons 
and the He Stromgren radius moves inward relative to r^. The function 
^gas^ other hand, depends strongly on both 6 and oi which, 

unfortunately, reduces its reliability in determining observa- 

tionally (see section 3 . 8 ). 

There are two classes of ions whose average ionization depends 

on (0 . Ihis dependence is defined arbitrarily as a variation of 
+ P 

(X ) greater that 10 % between models A and F. Class I 10 ns 
/+ + + + + +. 

(0, C,N, S,Ne and Ar ) all occur only at the Stromgren radius 
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where the variation in the spectrum of J ( t) is the greatest from 

V 

models A to F. As we discussed in the last section, the effect of 
albedo is to strongly enhance the intensity of UV photons near 
Hence the effect of increasing albedo is to enhance the lines from 
higher stages of ionization at the expense of the lines from lower 
stages of ionization. 

+ 2 +3 +2 +2 

Class II ions (C , C , S , Ar ) occur in ionization structures 

+ 2 

that are not dominated by the X ion. Dominance of the second ion 

+ 2 + P 

IS defined roughly as (X ) ^ 0.8. Physically changes in (X ) can 
+ PH 

only affect (X ) if the P + 1 ion does not dominate. The 
ionization structure of carbon is not dominated by any ion so all 
stages of ionization are affected by All of the effects are, 

however, much smaller than those of the variation of 5 . 

Although our central source temperature is for an HII region and 
not planetary nebulae, we can comment on the CIV/CIII ratio that has 
been a problem for models of NGC 'J02'J in particular (Bohlin, Marionni, 
and Stecher 1975 Pequignot, Aldrovandi, and Stasinska 1977) • 
Depending on the details of model, the model CIV/CIII ratio 
[ 1 ( 15 49 S/I ( 1909 ^) ] has been from 5 'to 10 times larger than the 
observed ratio. Our calculations show that increasing the albedo 
enhances the abundance of CIV relative to CIII. However, increasing 
6 has the opposite effect. It would appear that some combination of 
positive 6 and resonance scattering of the permitted CIV X 15^9 
line to increase its absorption by dust could help solve this problem. 

In Figure 6 we have plotted the ionization structures of H, He, 0, 
C, N and S for models A through F. The Ne and Ar ionization structures 
have not been plotted. The Ne structure is qualitatively like that of 
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SULFUR NITROGEN CARBON OXYGEN 



y/x 


Fxg. 6a. Tonizatxon Structures of He;, 0, C, N, and S. 

Model A xs the standard models Model B has § = 1 and 
Model C has 6 = 2. 
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D E F 



Fig. 6b, Ionization Structures of H, He, 0, C, N, and S, 

Model D has no dust. Model E has cu =0.5 and 

Model F has co = 1.0. 

V 
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O. However^ the Ar structure is unique as can be seen in table IV. 

+ 3 

It is almost dominated by (Ar ) and there is a monatomic decrease 
/ +P 

in (Ar ) as P decreases for the 6=0 models (C also exhibits 

this behavior) , The ionization structures are plotted versus t/t^ 

and there is no material with t > (see s.ec.tion 2.6). 

The ionization potentials of and are below l3,6eV so 

that they are not neutral at the H Stromgren radius, does not exist 
, 0 t 

in pur models but C and C co-exist at For clarity in 

^ 0 + 

Figure Oj C and C abudances are extended as straight lines for 
T > C ionization structure is the only one we have plotted 

that IS sensitive to the electron density. The discontinuities in the 
+ 2 

C structure in models B and C occur at the He Stromgren radius 

Sarazin (197?) Plo'fcs on his Figure 1 the ionization structures 
of He, N, O and S for a Salpeter initial-mass function central 
source whose spectrum is similar to an 06 star (T .s. 42,000°K, ZAMS) . 
This spectrum is slightly hotter and broader in the UV than our stan- 
dard model central source, Sarazin finds a monotonic increase in the 

* 

He transiton zone widths with increasing— 6 . Our results indicate 

that this behavior is not monotonic. This result is, however, highly 

dependent on the details of the central source spectrum and, to some 

extent, the method of solution. We find that albedo has little effect 

on the width of the transition zones, Sarazin 's zone II ions 

2 , "^3 ^3 "^ 2 . *^3 

(N ,0 ,S ,C ,Ne ,Ar ) are all dominant and have the largest 
average ionization in model A, The average ionization of these ions 
decrease with increasing §. The average ionizations increase slightly, 
however, with increasing albedo. 


-4, 

Sarazin uses a slightly smaller value of e (e ~ 1,4 xlO ) iov 
his models in Figure 1 which accounts, in part, for his wider 
transition zones , 
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Balick (I 975 ) calculated average ronizatxons using the on- 

*0 3 ~3 

the-spot approxxmation, T = 40,000 K, and n = lO'^cm for 

unxforra nebula and varxous forms for f ,(v) xncludxng a Ixnear form 

d 

+ P 

(6=1), He then varxes e and studxes the behavxor of (X ). He 

> 

fxnds that only an extreme step form of ^^(v) 

Vq s V < 1.807 vq, ^^(v) = 10 1,807vq «5 V g <x> ] can reproduce 

Mezger's (Mezger et al, 197^) observatxons of L. ^ ~ 0.1 for several 

He 

galactxc HII regxons. Our parametrxc studxes have shown that we need 
3 ^ o: ^ 4 to produce 0,1 wxth our standard model. 

Balick and Sneden (1977) have calculated model HII regxons using 
the on-the-ppot approximation with up-to-date model atmosphere central 
source spectrums that include metal opacities. ^ These atmospheres have 
considerably less high energy UV flux than those of Auer and Mihalas 
( 1972 ) or black bodies. Therefore, emitted low excitation lines will 
be enhanced relative to high excitation lines if more realistic 
stellar spectra are used. 

7 . 2.3 Infrared Emission 

In order to study the effects of albedo on the dust temperature 
structure, we will use the simplifying assumptions of generalized 
Mie theory (with j = 1) and optically thin IR emission. With these 
approximations the IR spectrum is a probe of the entire dust tempera- 
ture structure. We will use type II models so that the luminosity 
of the central source is constant. 

Figure 7 shows the IR spectrum for our standard model with 

oj =0 and a> =0.9 independent of frequency. The spectrum was 
V V 

calculated using Eq. (4.21). The albedo therefore applies only to 
the UV and OP photons that heat the dust. Again we assume that there 
IS no material outside of the Stromgren radius of the gas . If 
radiative transfer of IR photons had been calculated, then because 
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SlR(iJ) (lO’^®ergs/cm^-sec-HE) 




most of the spectrum would be optxcally thin due to the X 

dependence of K (r) in the infrared, the main effects show up in the 

near infrared where the unattenuated Rayleigh-Jeans tail of the central 

source adds significantly to S^j^(v) . 

There are three important effects illustrated by Figure 'f . First, 

as the albedo increases the luminosity of the IR spectrum decreases 

(i.e.lC^(r) is smaller for a given extinction coefficient so that the dust 
V 

absorbs less energy) . Second, the spectrum is shifted towards longer wave- 
lengths and/ IS reduced in width (i.e. the spectrum looks more like that of 
a diluted black body) As we have seen, dust scattering transports photons 
outward and thus tends to make the mean intensity more uniform as a function 
of radius (see. Fig.4) . Hence the effects of albedo are 1) to reduce the 


maximum dust temperature which occurs at the center of the nebula 


(i.e. T, 182 K for a> = 0 and = Il7 K for co = 0*9) and, 

d,max V d,max ' ^ 


2) to reduce the variation in dust temperature throughout the nebula 

(i.e. T, /T, = 4.75 for OD = 0 and T /T = 4.02 for 

CD = 0 . 9 ). The first of these effects is due primarily to the (1 - cd ) 

V V 

term in the right-hand side of Eq. (4.4) while the second of these 
effects IS due primarily to radiative transfer. Thus we see that as 
albedo increases, the dust temperature at the same relative radius 
decreases and, as a result, S (\>) for the near and raid infrared is 
substantially lower than the zero albedo model. Third, because <r^ 

IS larger for the non-zero model, there is more material with low dust 
temperatures which radiates the "excess" far IR flux. 

7.3 Comparison of Solutions of Radiative Transfer 


In this section we will make detailed comparisons of the solutions 
of radiative transfer that were discussed in chapter 5 with the QDM 
solutions and with the Eddington approximation solutions. We will 


concentrate on observable parameters such as f ^ ^ » the average 

ionizations and the IR spectrum. However, the hydrogen ionization 
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structure will be examined to extract the differences in the properties 
of the approximations. 

7 . 3.1 Diffusion Approximation 

On Figure 8a we plot the emergent luminosity normalized to the 
luminosity of the central source for pure dust nebulae with isotropic 
scattering as a function of the total extinction optical depth 
(f^(v)=l). The diffusion approximation (see section 5*3) does a very 
good job of predicting the emergent luminosity when the albedo is 
large and ^ 2. The free surface diffusion approximation is almost 
exact for co^ ^ O. 9 , whereas the infinite diffusion approximation is 
10 to low even for cOq ^ 

The diffusion approximation is not valid at the boundaries and 
the free surface boundary condition uses a calculation similar to QKW 
to correct the emergent luminosity. We therefore expect good agreement 
when the other assumptions of the diffusion approximation are met 
(see section 5.3). A more exacting test of the accuracy of the 
approximation is a comparison of L(t)/L* when all assumptions are met. 
Figure 8b plots the transmitted luminosity as a function of t for 

= 10 and CD = O. 9 , O. 95 , 0.99^ and 1.0. 
f V 

The QDM solution for cd = 1.0 is a test of the overall accuracy 

V 

of the differential equation solvers of program SCATER. We mentioned 
in section 6.5 that flux was not strictly conserved in this case. 

Over most of the volume of the nebula, the flux is 1.5?§ high, a result 
we attribute to the linear approximations made in doing angular 
integrals. This error increases to 3^ at the center of the nebula 




88 




Fig. 8a. Emergent Luminosity Versus Tt For Various o> and For the 
Diffusion Approximation and the QDM. The solid lines are 
for infinite nebulae and the dashed lines are for nebulae 
with a free surface. 
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Lm/L 



T/'f, 


Fxg. 8b. Comparison of Luminosity Versus t for the Diffusion 

Approximation and for the QDM. = 10 for all models. 
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where the angular grad is coarse (i.e. there are only 10 grid points 
at r =: whereas there are 110 at r = r^ if all of the radial arrays 

are used) and where the r dependence of ( t) results in less 
accurate radial integration. 

As CD decreases, the diffusion approximation becomes more accurate 
V 

over most of the volume of the nebula. For cd < 0 . 99 } the infinite 

V 

and free surface approximations differ only at the outer boundary, and 

agree with the QDM solution to within 10^ over most of the volume. 

When CD = 0.99^ the infinite and free surface solutions differ by 20^ 

V 

while the QEM and free surface solutions agree to within 0.5^. Part 
of the large differences between the QDM and diffusion solutions at 
the center of the nebula must be attributed to the inaccuracies of 
program SCATER. 

7.3.2 Average lopizations of Type I Models 

Ihe accuracy of any approximate solution of radiative transfer 
must be assessed in terras of observable quantities. To do this, we 
will define the QEM solution to be an "exact" solution. We will cal- 
culate models A and F (see section 7.2.2) using the Eddington approxi- 
mation, the modified Eddington approximation, and GOSA (model A only) , 
Table V lists the observable parameters for these models and for the 
various solutions . 

For model A average ionizations, both the Eddington approximation 
and GOSA differ from QIM by less than 1.5?^ all ions, and f 

gas 

vanes by less than 2^ among the solutions. For model F however, the 
average ionizations can vary by as much as ^5-5^ in one case (Ar^ ) . For 
dominant ions (see section 7*2.2), the Eddington approximation values 
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TABLE V 


IONIZATION STRUCTURE OF MODEL NEBULAE 
COMPARISON OF SOLUTIONS OF RADIATIVE TRANSFER 

■H" 

Models 


parameter 

A 

A 

A 

F 

F 

F 

Solution 

edd’’ 

GOSA 

QDM* 

edd"^ 

QDM* 

MEA^ 

^0 

4.128 

4.128 

4.128 

1.300 

1.300 

1.300 

6 

0 

0 

0 

0 

0 

0 

CD 

\\ 

0 

0 

0 

1.0 

1.0 

1.0 

V 

f 

gas 

0.298 

0.291 

0.298 

0.961 

0.936 

0.961 

^He 

1.015 

1.016 

1.016 

1.014 

1.016 

1.015 

<x) 

0.973 

0.972 

0.973 

0.976 

0.970 

0.973 

(y> 

0.988 

0.988 

0,988 

0.989 

0.986 

0.988 

<o'"> 

0.035 

0.034 

0.035 

0,027 

0.032 

0.030 

<o"^> 

0.952 

0.953 

0.952 

0.962 

0.953 

0.957 


0.043 

0.o43 

0.044 

0.036 

0.046 

0.042 

<c^"> 

0.426 

0.420 

0.425 

0,375 

0.385 

0.375 

<c-"3) 

0.525 

0.531 

0.525 

0.586 

0.564 

0.580 

<N^ > 

o.o43 

0.042 

0.043 

0.035 

0.043 

o.o 4 o 


0.818 

0.816 

0.817 

0.819 

0.809 

0.801 

<N^^> 

0.135 

0.138 

0.135 

0.142 

0.142 

0.154 

(s'") 

0.0097 

0.0098 

0 .0098 

0.0083 

0.012 

0.0097 


0.068 

0.066 

0.068 

0,051 

0.062 

0.058 

<s-^ 3 > 

0.733 

0.732 

0.733 

0.735 

0.723 

0.715 


0.188 

0.192 

0.189 

0.205 

0.202 

0.217 

<Ne' > 

0.133 

0.130 

0.133 

0.104 

0.121 

0.115 

<Ne"^> 

0.864 

0.867 

0.864 

0.893 

0.876 

0.882 
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Models*"'' 


Parameter 

A 

A 

A 

F 

F 

F 

Solution 

EDD^ 

* 

GOSA 

QEM* 

edd" 


MEA=^^ 

<Ar > 

0.0056 

0.0055 

0.0056 

0.0046 

0.0064 

0.0054 

<a /5 

0.29^1 

0.288 

0.294 

0.243 

0.266 

0.257 

<a/ 3 > 

0.698 
1 

0.704 

0.698 

0.750 

0.725 

0.736 


GOSA; Generalized On-the-Spot Approximation (Section 5 . 2 . 1 ) 

EDD- Eddington Appro'ximation 
^DM* Quasi-Diffusion Method 

1 %EA- Modified Eddington Approximation (Section 5 -^) 

Models ^EDD reproduced from Table IV for comparison. 


93 

















agree wxth the QEM values to withxn lO^ xn all cases. In general, 

the modxfxed Eddxngton approxxmatxon (mEA) does better than the 

Eddxngton approxxmatxon, especxally for those xons whose average 

+2 

xonxzations are small. Only for (C ), (n and (S ) does the MEA 

faxl to do a better job of reproducxng the average xonxzatxons. 

7 . 3.3 lonizatxon Structure of Type II Models 

The detailed lonizatxon structure of hydrogen xs presented xn 

Exgure 9 £ 0 ^ type II models. The solxd Ixnes are for QEM models with 

= 0, 0 5^ 0.8 and 1.0 xn the UV. In addxtion we plot (1 - x)* for 

the GOSA, the Eddxngton approxxmatxon and the analytxc model. For 

sxmplxcxty helxum has not been xncluded xn these calculations although 

all the other parameters are the same as our standard model. 

For 03 — 0 we can make a direct comparison of all our models. As 

V 

can be seen from Figure (1 “ ^) calculated wxth GOSA and Eddington 
approximation agree with QDM's calculation over most of the nebula. 

The Eddxngton approximation agrees to within 5 ^ over all of the nebula 
and to within 2^ except where (1 - x) /... 0.1,. GOSA agrees to within 
105 ^ when (1 - x) 9^ 0.1 but varies by as much as hO^ when (1 - x) Rs 1.0 
All three models produce exactly the same value of (x.e. they all 

find the Stromgren radius condition [Eq. 3.32] satisfied at the same 
radial grid point). Ihe analytic model which uses OSA (diffuse UV" 
photons are absorbed on-the-spot by the gas) underestimates the value 
of (1 - x) because dust does absorb a large fraction of diffuse UV 
photons. The analytic model also underestimates the value of 
by 10^ primarily because it uses frequency averaged intensities and 
cross sections. Hence the analytic model cannot account for the fact 
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NEUTRAL FRACTION OF H 



Fig. 9 . Neutral Fraction of Hydrogen Versus t for Various 
and Solutions of Radiative Transfer. A model 
without dust ( dashed line) calculated with the QDM 
IS presented for comparison with the = 1 model. 
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that hard UV photons penetrate farther xnto the gas than do photons 
with V ^ Vq- rounding of 1 - x at the Stromgren radius xs due to 

ionizations by recombination lines and by penetrating hard UV photons. 

As the albedo increases the Eddington approximation becomes less 

accurate when compared to the QDM solution of (1 - x) because the 

Eddington approximation does not account for the forward peaking of 

the scattered continuum. Ihe Eddington approximation model calculates 

a larger value of J^(t) over most of the nebula and hence a smaller 

V 

value of (1 - x) . It should be remembered from the last section that 
despite these differences in the internal structure^ observable 
quantities agree well between the two models , 

We have also plotted 1 - x for a dusty model with oij = 1 and a 

V 

dustless model [dashed line with dq- = (2 ^ 10 )n (r)crQdr] with the 
same luminosity and gas density. The difference in (1 - x) between 
the two models is due primarily to the difference in the stellar flux 
as a function of radius . The dusty model quickly scatters the stellar 
photons into diffuse (and more isotropic) photons, increasing the mean 
intensity near the center of the nebula and thus reducing (1 - x) 
relative to the dustless model. At the StrJ5mgren radius the dustless 
model has a smaller (1 - x) and therefore a larger because the 

stellar flux is considerably larger than at a corresponding point in 
the dusty model. Hence in the dustless model the stellar photons 


penetrate farther into the nebula before they are absorbed by the gas. 

7.3.^ Variation of f with Albedo 
gas 


The fraction of UV photons "destroyed" by gas, depends 

primarily on the value of for our standard models because of the 
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"sharpness" of the ionization fronts. We therefore expect that all 

models of radiative transfer that accurately calculate will also 

accurately calculate ^ gg^g • On Figure 10 we plot the variation of 

fgas (ca-lculated with Eq. 3.3^) with albedo for the QDM,- Eddington 

approximation, GOSA and analytic models. As expected the Eddington 

approximation calculation of agrees very well with the QDM 

calculation. In addition, both GOSA and analytic models accurately 

calculate f when o: =0. 

gas V 

Although we have not plotted L. , the ratio of the intensity 

He 

of helium to hydrogen recombination line (Eq. 3.3^), we find similarly 
excellent agreement among the various models because this quantity 
also depends primarily on the value of 
7.3-5 Infrared Emission 

On Figure 7 we have plotted the infrared spectrum for our Standard 

models and for the QDM solution, the Eddington approximation, GOSA, 

and the analytic model. For to = 0 all approximate solutions reproduce 

V 

the shape and the strength of the spectrum. 

The analytic model is least accurate for two reasons. First, 

the volume integral of Eq. (2.21) is calculated ma more approximate 

way in the analytic model. Second, the analytic model underestimates 

Ti l>y and hence underestimates the far infrared flux that is 

radiated by low temperature dust at the edge of the nebula. 

Agreement between the Eddington approximation, GOSA and the QDM 

spectra for oi = 0 is almost exact because, for this model, most dust 

heating is by stellar photons in the OP and UV ranges. However, for 

CO > 0^ the Eddington approximation underestimates at the center 

V V 
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FRACTION ABSORBED BY GAS 



Fiff. 10. Variation of f Versus m for Various Solutions 

gas 

of Radiative Transfer. 
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of the nebula and overestimates at the Stromgren radius. Hence 

V 

the near infrared is too low compared to the QDM and the far infrared 

IS too high. Even so, L (Eq. 4.23) agrees to within 0.2^ 

I R 

between the Eddington and QDM spectra. 

In general, these comparisons are academic because of the 20 to 

50^ errors in the published IR broadband photometric data. 

7 . 3.6 Accuracy of Semi -Analytic Solution 

The semi-analytic solutions described briefly in section 5*2.2 

and appendix B and completely in Papers I and II are extremely useful 

for calculating f and D , because they are simple to calculate. 

gas ^ + 

The models use frequency averaged intensities and fluxes, however, so 
that they cannot be expected to accurately describe the details of the 
ionization structures at the Stromgren radius. 


We have run several GOSA calculations using program SCATER with 


various power law forms for the frequency dependence of the dust 
extinction coefficient in the UV ranging from (v/vq) to (\>/vq) > 

with 1*2 the range 1 ^ 3^ with Yq == 0.04 and O.IO 8 (see 

Eq. B.l), and with Y = 0.1. We have also calculated the semi- 

analytic models with values of O (see Eq. B.4) calculated by using^ 
the corresponding power law from for c^(v) averaged over the spectrum 
of the central source [i.e. in Eq. B.2, replace o'jj(v) with cr^(v) 
and use J^(r^) ] . 

When we compare these solutions, we find that for (X ^ 0 


(6^0 in Eq. 7.2) 
within 20^ in all cases, 
agree to within ^0^ but 


L agrees to within 10^ and f to 
He’’" 

However, when o:>0(8>0) L , usually 

He 

f can differ by as much as 150^. The 
gas 
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reason for these large errors is that the semi -analytic models fail 


to predict the sharp decrease in as a positive 8 (defined 

ahove') increases. Physically what is happening is that because the 

semi-analytic model treats all tJV photons alike^ it fails to account 

for the ’’penetrating power" of hard UV photons (remember that the 

gas opacity drops rapidly as v increases) . If these photons are 

preferentially absorbed by dust, then the Stromgren radius must move 

inward. Even so, L is fairly accurate because it is calculated 

He 

as the ratio of two functions of (see Eq. B.33). ^gas^ 

other hand, is proportional to a single function of (EQ* B.31). 

Finally, we can define an effective optical depth, which 

when used in Eq. (B.31)> produces the same value of f that is 

gas 

calculated by program SCATER. If we do this as a function of the UV 
albedo for our standard model parameters, we get the curves shown in 
Figure 11. Here we have plotted t ~„/t' where t — (1-o>)t-, 
versus co. We see that these curves bear little resemblance to those 
of Figure 7 paper I . If, however, one were to continue the lengthy 


and costly process of calculating type I models (S^ must be adjusted 
to keep t 2 _ constant as co increases), then Figure 11 could be used 
to adjust the values of deduced from observations (section 3,8) 

for an assumed albedo. 

7.^ Conclusions and Concluding Remarks 

We have considered the effects of dust scattering on the internal 
structure and observable quantities of dusty nebulae. We have also 
compared our solutions with more approximate solutions of radiative 
transfer. When doing these comparison we consider the QDM solution 
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to be exact Particular emphasxs has been given to ionization 
structures of the gas and to the emitted infrared spectrum. The major 
new results and conclusions of this study u.re listed below, 

1. UV and OP recomhination line intensities ( even those not reso- 
nantly scattered by the gas) remain fairly isotropic throughout 
the nebula while the scattered continuum intensities become 
outwardly peaked. 

2. The radial variation of the Eddington factor is fairly indepen- 
dent of albedo. 

3 . does not depend on the albedo of dust. It does depend 

weakly on and strongly on the frequency dependence of 'the 

dust extinction coefficient and on the spectrum of the central 
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d) In all cases the effects of the variation of the albedo 

were much smaller than the effects of the variation of the 
frequency dependence of the dust UV extinction coefficient. 

6. Increasing the albedo of dust in the frequency ranges primarily 
responsible for dust heating has the following two major effects. 

a) The total infrared luminosity and the dust temperature at 
the same relative radius are reduced because the dust 
absorption coefficient is reduced for a given extinction 
coefficient . 

b) The variation of the dust temperature from the inner to the 
outer radius of the nebula is reduced, thus making the 
observed IR spectrum narrower. 

7. The diffusion approximation with the "free surface" boundary 
condition can accurately predict both the emergent luminosity 
and the internal variation of the flux when the assumptions of 
the approximation (large scattering optical depth) are met. 

8 . When the UV albedo is zero the average ionizations of all elements 
calculated using GOSA or the Eddington approximation agree with 
the values calculated by QDM to within 1.5??. f vanes by less 
than 2^ among the different solutions . 

9. When o> =1 for UV photons the Eddington approximation average 

V 

ionizations agree with the QDM values to within 10^ for all 
dominant 10ns . Variations between the solutions as large as 
can be found when the average ionizations are small. In 
general, the modified Eddington approximation predicts values 
closer to those of QDM. 
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10 . 

11 . 


^gas specti‘um are accurately calculated by the 

Eddington approximation . 

The semi-analytic solutions calculate L . and f to 

He* 

Within 1055 and 20 ^ respectively if dust does not preferentially 

absorb hard UV photons . when this does happen^ L ^ is 

He 

accurate to within 50^ but f can be off my more than a 

gas ^ 

factor of 2 . 


Finally, a few words should be said about the conclusions we did 
not or could not make,, As has often been the case of previous genera- 
tions of graduate students, this thesis is being completed at this 
time because a new and completely different career awaits in the very 
near future. Hence several things must be left undone. 

Although program SCATER has been used in close collaboration with 
the author by Alan Tokunaga of NASA-Ames Research Center in an attempt 
to fit the IR spectra of SlUo, this work is still too preliminary to 
include here. Regretfully, we have not had time to use SCATER to fit 
both the OP and IR spectra of an HI I region. 

In addition, we have not had time to include radiative transfer 
of OP and IR forbidden lines . This work would require first a self- 
consistent electron temperature calculation and then calculation 
of the populations of the collisionally excited states from which 


these forbidden lines coiie. 



Appendix A 


DIFFERENCE EQUATIONS 


The JAY [Eq. (2.19)] and RAY [Eq. (2.^2)] equations are both of the 

form 


where the known functions A through F are, in general, dependent on 
frequency, x , and Y . The boundary conditions are of the form 




= (HY + l) 


x=x 


(A.2) 


x=x_ 


B 


for both the inner and outer boundaries of the JAY and RAY equations. 
The range of x is divided into n grid points with 

^1 ^2 ^ *^^i 1 ^ 1+1 ^n* ^ three-point 

central difference approximation for the jgY term evaluated at the 
th 

1 grid point gives 


B ^ 
1+1 + 

B ' 
1 


D ^t-D y1 

1+1 1+1 1 1 

* B 

1 



’D Y -D tY ,1 
1 1 1-1 1-1 

C 

1+1 

C 

l_ 


i 1+1 1 J 

_ 1 

^-xj 


~X - X 

L 1 i-1 . 


1+1 ' \-l 


(A. 3) 


Substituting Eq (A. 3) into Eq (A.l) gives the set of difference 
equations that approximate Eq. (A 1) • 
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( V V V “i"- 

+ (’'3^ V - Vi 


(A. 4 ) 


where 


T, = 


C C B , 

A 

1 

1,1 

1+1 1 1-1 

1 


w. 

c c , B 

A 

2 

1,1 

1+1 1-1 1 

1 

= 


C C , B , 

A 

3 

2,1 

1 1-1 1+1 

1 

II 

c 

2,1 

' - C T ® A 

1+i 1-1 1 1 


C , ^ C C ^ 

5 3^3. 1+1 1 1-1 


(A.5) 


for 1=2, 3 , 
defined as 


W. 


n- 1 . The weighting functions W. are 




W, 


2,x 


1+1 1 


X - X T 
1 1-1 


(A.6) 


w 


5,=(x,,-x) (x,T-x ,) (x-x ,) 
3,1 1+1 1 1+1 1-1^ ' 1 i-l'^ 


Note that the W functions depend only on the grid spacing and hence 
can be calculated when the grid is defined and/or changed and stored 
for later use. Using the Eq. ( 2 . 33 ) and Eq. ( 2 .h^) forms of the 
JAY and RAY equations respectively, the W functions depend only on 
T(r) and can be used for both equations. For the RAY equation p is 
tabularized as a function R and B. 

At the boundaries DY is expanded in a second-order Taylor series 
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in a manner described by Auer (I 967 ) . Eq. (A. 2) is used to evaluate 
the first-order derivative of DY and Eq. (A.l) is used to evaluate 
the second-order derivative. At the inner boundary we get the 
following relationship between and Y^ 




Y1H1 + 




W2^2^1^2 ^ - 




1 ^ w 


(A.T) 


(inner B.C.) 


where 


h - Kv ■ 


(A.8) 


At the outer boundary. 


-ABC tD Y ^ + 
n n n-1 n-1 n n-1 


[A B. 


^ C ,D G - 
n n n-1 n n 


V H + 

'n n 


6 F 1 
n n-' 


Y = - 
n 


6 F + 
n n 


^n^n 


(A.9) 


(Outer B.C.) 


where 


Y = 1a (X - 
’n n n 


X J (3B C T 
n-l"^ ^ n n-1 


- B tC ) 
n-1 


6 


n 



X , ) 

n-1 


2 


C C G 
n n-1 n 


(A. 10) 


The inner boundary of the RAY equation requires more care because |i 
is identically zero on the symmetry line. Thus we use Eqs. (2.^2) 
and ( 2 . 53 ) at the inner boundary. 
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Eqs. (a. 4 ), (a. 7) (a. 9) form a diagonally dominant tri- 

diagonal system of linear equations that are easily solved using 
the Gausian elimination method (Greenspan I97O) ♦ 

Finally, the W weights can be used to calculate radial 
derivatives of a general function f(x) using a three-point 
approximation* 


df(x) 

dx 






(A. 11 ) 
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APPENDIX B 


SEMI -ANALYTIC MODELS 


B.l Introduction 

In this appendix we discuss the ionization structure of hydrogen 

and helium. This discussion is a condensed form of Petrosian and 

Dana (1975^ Paper l) and Dana and Petrosian (197^^ Paper II) For a 

general form of the frequency dependence of the dust extinction 

coefficient in the UV, we derive coupled differential equations which 

can easily be solved numerically for the Stromgren radii of hydrogen 

and helium. We present approximate analytic expressions for the 

evaluation of the variation of the fraction of ionizing radiation 

absorbed by gas, f .and the ratio of the volume emission measures 

gas^ 

of Hell to HII regions, L } with the spectrum of the ionizing 

He*^ 

source, helium abundance and absorption properties of dust As before 
we assume spherically symmetric nebulae and we use frequency averaged 
cross sections and intensities. We will apply these models to dusty 
HII regions excited by 0 or B stars so that no doubly ionized 
helium exists. Equivalently, we could apply these models to the 
regions of the nebulae outside of the Hell I Stromgren sphere. 

B.2 Equations of Transfer for UV Photons 

We shall be dealing with fluxes and intensities integrated over 
frequency. We must distinguish between photons capable of ionizing 
only hydrogen and those capable of ionizing both hydrogen and helium 


1Q9 



We therefore define net fluxes through spherical shells 

fo 

S^.(-r) = / dV S^(r), S^{r) = / dv S^(r), 


"0 


1.807V 


0 


(B.l) 


S(r) = Sj_(r) + Sg(r), y{v) = S 2 (r)/S(r), Yq = 


and similar expressions for the intensities iCr) and J(r), We 
use 4 vq as the maximum frequency because photons with V > 4 Vq 
will be absorbed in a small Helll region interior to the region 
of validity of these models. 

Average hydrogen cross sections for V < l.SOYV^ and for 
V > 1.807 Vq are defined as 


1.807V, 


0 


/ dv ajj(v) J (r), 


'0 

4v 




1.807 v ^ 




(B.2) 


so that the total average hydrogen cross section is 




1+ 

1 + Jg/Ji 


(B.3) 


^ ~ ^ " %2^'^H,1 * 


1.10 



Similarly, for dust we define cr^ ^ and 2 that 


*^d “^d,! 


1+ (l+a)j2/J^ 

1 + Jg/J^ 


, 1 + « • ■ (B-4) 


In these models we will characterize the dust by its optical depth 
averaged over the frequency range Vq ^ V < 1.807vq* Hence the dust 
extinction cross section at the Lyman limit is replaced by a 

a, -L 

A 

in our formulas for (Eq. 3.24), e (Eq 3 28) and for the optical 
depth functions t and R. 

For helium we have one average cross section (the He 

notation is unnecessary because does not exist in this model) 

defined as 


4v, 




I.BoTv 


0 


(B.5) 


We now define separate average cross sections a and or for the 

X X! 

stellar and diffuse radiation. For the average intensities 

* * 

and can be replaced by and , respectively. 

We will use the on-the-spot approximation (i.e. diffuse UV 
recombination photons are absorbed by the gas where they are emitted) 
and set the dust albedo co =0. With these approximations the diffuse 


UV radiation satisfies the relations 



(B.6) 


(1 - x)n(r)cTj^^^ kit JD (r) = xn(r)n^ aj + g^Yyn(r)n^ 

(1 - y)Yn(r) o-jj^ (1 + (r) = Yyn(r)n^O!^g 

where xs the fraction of photons from recombinations to excited 

states of helium that are capable of ionizing hydrogen 
with the approximations of Section 3.^) and 


P = - x)ajj(1.807VQ)/[Y(l - (B T) 


IS the fraction of recombination photons to the ground state of 
helium that are absorbed by hydrogen . Upon substitution of 
equations (B.6) into the ionization equilibrium equation (Eq. 3 2 
with 8=6^=: 0), the dimensionless ionization equations (Eq. 3.29) 
become 


* 




(1 - x) = (3d-Yy)[x-yY(S^+C2)] V ^ ® 


(B.8) 


(1 - y) = y(xfYy)[l y( t) S* ( -r) ] 


In most OSA calculations this fraction (and consequently the quantity 
S in Eq. B.8) is set equal to zero. If this were the case, the 
absorption of photons by hydrogen could also be neglected. As 

we shall see below, neglecting p or ^2 with respect to unity will 
cause up to a 30 percent underestimation of Yq ^ similar over- 
estimation of Y. 
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* 

where 




V - =“He ^ 


(B.9) 


and we have defined 




(B.IO) 


Substitution of these equations into the equations governing radiative 

transfer of stellar photons ( dS /dr = -K , , S ) gives 

^ 1 tot, 11 ^ 


dS*/dr = -4rtr^n(r)n^0^ [x+ Yy(C£<- Cg) ]cfH*i 


* * * * 


dSg/dr = -'^^^ 2 . ®2 I ^2^^^ '*’ 


* -jf- * • 3 '- 


[X - yy(Sj_H- C^)] <^h, 2V‘^H^ j 


1 . 


(B.12) 


Addition of these two equations gives 


•X- 

dS 


/dr -IC*S*-l^rtr^n(r)n^O^ [x + yY(l - Ci)] , (B-13) 


1 ' 


which for =1 is identical to Eq. ( 25 ) of paper I for nebulae with 

* -X- * 

dust and hydrogen only (note also that K^/K^ ^ = 1 + CC y depends on y) • 


*Note that cr* = <^h*1 
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(From here >on we shall be dealing only with stellar photons. 
Therefore we shall elitninate the asterrsk notation in what follows. 
If we eliminate and m favor of S and y in lEqs 

and (B.12 )j Tthe general equations for stellar , photons .are 


A 

ds/dT = -(1 + ay)s - 2(t) 


■dy/dT = -0!Y(1 - y) + (1 - y')Oy-^ )(i!- - Py)“^ 2(t)/S 




.where 


dT = K^^^dr, S(t) = S (t)/S^ 

.i:(!r) = [x+ ty(i - C;i)j(-x-+ yt)![r^g|/g^]/Sq (B. 15 ) 

yY(l + iQj 

y' = 

X + yY(i - 

'We shall Ignore •(! - xi) and (1 - y) with respect to unity m 
Eqs. (lB.il) - (B.I 5 ) (te.'g. -paper j) . This is true 'everywhere .(m 
particular for large dust optical -depths) -except wery near the 
ronisation fron.ts,. This amounts to setting u = i ,and y = 1 
inside of their respective Stromgren spheres in Eqs. {3,. 11) — ’,(©. 15 ). 



B.3 Pure Hydrogen Kebulae 


For nebulae without helium (y=Y = 0:=p = 0), Eqs. (b 8 ) and 
(B.14) reduce to 


(1 - x) = x^A^R^Gg/S(T) 

(B 16) 

dS/dr = -S - x^(R^Gg/G^)/SQ 

(B 17 ) 


These equations are equivalent to Eqs. (23) and ( 25 ) of paper 
Eq (B.I 7 ) IS easily solved to give 


S(<r) = 1 - g(<r)/SQ] (b.18) 


where 


g(T) 


j dt xV(t)Gg(t)eVG^(t) 
0 


(B.I 9 ) 


and we have used the normalization of S(t) (Eqn. 2.30) s-s the initial 
condition of Eq. (B.I 7 ) . In practice x is set equal to 1 inside 
of the hydrogen Str'dmgren sphere and gC^) becomes a simple function 
of <f once the radial dependence of the dust and gas densities 
[Gjj(r) and G^(n)] have been specified. For uniform nebulae. 


g(t) = g(t+ Rq) - gCR^) 



2 


X 


2 x+ 2 


(B.20) 
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where is the inner radius of the nebula in optical depth units 

(Eqn. 2.24). 

A 

Equation (B.18} points out the role of Sq in determining the 
optical depth of a nebula. The Stromgren radius is defined in this 

A A 

model when S(t) = 0. This occurs when g(T 2 ^)= S^. (Note that g(T) 
is a monotonically increasing function of t for the power law 
density variations of section 2.-4) . Although this is the only model 

A 

where can be detennined analytically^ plays an analogous 

role in all of our models including the QDM model. 

The ionization of hydrogen can now be calculated from Eqs. (B.16) 
and (B.18), Interior to the Stromgren radius, (1 - x) is given by 
the right-hand side (with x = 1) of Eq. (b. 16). At the Stromgren 
sphere, however, the quadratic form of Eq. (B.16) must be solved to 
give 


(1 - x) = [1+ A(t)1 - V[1+-A(t)]^ - 1 


(B.21) 


where 


A(t) = S(t)/2A^V . (B.22) 

o 


For nebulae without dust, Eq. (B.I 9 ) with ~ 0, gives 

dS/dTjj = , dTjj = (B. 23 ) 

A 

and Sq is given by Eq. (3.24) rath hjj(3^o)cr^ q = X* The 

solution to Eq. (B.23) is 
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(b.24) 


S(Tj,) = 1 - ?(Tg)/So 


where 

g(TH) = / dtxV(t)G^(t) (B.25) 

and <r„ 2 .S the optical depth of neutral hydrogen at the Lyman 

H 

A 

limit. Again, the role of is apparent from Eq. (B.24). For 
uniform nebulae (with x = 1 inside of Stromgren sphere), 


?(t„) - iWjt 



(B.26) 


The ionization structure is given by Eq. (B.21). 

B .4 Nebulae with Hydrogen and Helium 

Since to our knowledge there are n o approximate solutions for 
nebulae with hydrogen and helium but without dust, we consider this 
case first. With = 0, Eqs. (b. 14) and (B.I 5 ) reduce to 


S(Tjj) = 1 - i('Tjj)/SQ , 

(B.2T) 

^ _ zi (1 - v)(Pv - yQ 

A A 

dS S (1 - Py) 
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where 





dt[x+ yY(l - S^)](x+ Yy)R (t)Gg(t) 


(B.28) 


0 

To evaluate Eq (b.28)^ x and y are set equal to 1 xnside of their 
respective Stromgren spheres and are set equal to 0 outside. 

The parameters p and Y^ in Eq (B,28) vary throughout the 
nebula The variation of p is due to the change in the spectrum of 
UV photons. However j since the photoionization cross section of 
hydrogen decreases rapidly with frequency (CTjjCy) « v 
P « 1 - ( 1 / 1 . 807 )^ ^0.8 and changes by a few percent for a variety 
of plausible spectra (cf . Table VI ) We therefore neglect the 
variation of p. The parameter y' ; on the other hand, varies because 
of the variation of which is primarily due to the change of the 

ratio (1 - x)/(l - y) in Eq (B. 7 ) . In general, if Yq < 

« 1 and it can be ignored However, for Yq > C 2 ^ 
the inner regions and near the edge of the n ebula where y 

-♦ /o^ ^ 0.6 (for T = 10^ °k) . As we shall see below, even 

in this case the variation of is negligible. Thus we shall assume 

that Y' is also constant. 

Equation (B. 27 ) is now readily solved* 


S(Th) 


/I-y(Tjj) 

\ 1-Yo 



1 -Y^ 

/ Pyq--^ 


(B.29) 


Figure 12 shows the variation of y with S for various values of 
Yq (solid lines) For Py^ = y' , y = Yq = constant. For Py^ < Y" , 

A 

y (and therefore S^) becomes zero at 
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TABLE VI 


AVERAGE CROSS SECTIONS AND STELLAR PARAMETERS 


T( °K) 

Yo 

6 

^H,1^°'0 

'^He'^'^O 

90300. 

0.500 

0.837 

0.452 

0.688 

62200. 

0.300 

0.808 

0.487 

0.832 

50000. 

0.194 

0.793 

0.514 

0.900 

1^0000. 

0.108 

0.787 

0.555 

0.959 

40000- ++ 

0.268 

0.831 

0.528 

0.802 

37500. 

0.0887 

0.786 

0.564 

0.979 

30900. 

0.0445 

0.785 

0.602 

1.023 

30000. tt 

0.0013 

0 

0 

0.747 

1.123 

n= 2* 

0.553 

0.916 

0.507 

0.410 

n=4* 

0.169 

0.855 

0.585 

0.754 


Averaged over power law distrxbutioix S =S^(v^/v)”. 

V U 0 

Averaged over black body distributxon B (t). 

V 

^^Averaged over model atmosphere (Auer and Mxhalas 1972) 
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Fxg, 12. Variation of y(t) Versus S('f) for Various Tt and y . 

Hie values y =0.^0} O.X^h and 0.0445 correspond to 

black body stellar temperatures of 90 ^ 300 } 5^^000 

and 30^900 °K respectively. The solid lines are for 

nebulae without dust (t 2=0). Four different dust cross 

sections were used* Ofe -1 (dash-dot line), CC = 0 (short 

dash) , Ct = 1 ( dot) , o; = 5 ( long dash) , where 1 + a = 

O', n* Two cases where t-,= 2 are plotted and labeled. 

^ Q.^ J. ^ 

All other cases are for t^= 1. All curves are for uniform 
gas and dust distributions^. 
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A 


S = 


A 


s 

cr 


(I-PYoA- ) 


(B.30) 


or at r = where r^ xs obtaxned from Eq. (b 28) with 

A 

I = s 1 - . Thus in these cases the helium Stromgren sphere 

IS inside of the hydrogen Stromgren sphere. For ^ 

Y increases toward the outer edge and approaches unity at the hydro- 

A 

gen Stromgren radius where S = 0. As is evident; the shape of these 
curves is determined primarily by the value of Yq (actually by the 
value of ^ * 

A 

Once the variations of S and y known; the ionization 

structure can be calculated with Eqs (B 8 ) - (B.10), Since 

2 

^ < (C^ ^ (B.8) can be approximated as (1-x) = A(t)x 

(A(t) is equal to A times the radial functions on the right-hand side 

H 

of Eq. B.8 ) a few values of a and a /c^ are given in 

-U O Xx0 O 

Table VI . The slow variation of these quantities due to changes in 
the spectrum of TJV radiation throughout the nebula are neglected in 
this treatment. 

Equation (B. 8 ) can now be solved for (1-y) using the above values 
of X and assuming C 2 “ constant. This assumption is clearly justified 


because R 

variation of results of these calculations are shown in 

Figure 13a for various values of Yq^ ^ = 0.1; and uniform nebulae. 

We also show the variation of As is evident; the assumption of 

constancy of is a good approximation for small values of y^. 

For small values of Yq^ ^2 ^ that y' Y ni I.O 5 Y (for 

T = 10^ °K, c.f. Burgess and Seaton i 960 ). For large values of 
0 


A 

(1 -Py)/S varies much more rapidly than any expected 
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'" o.o 0.2 0.4 0.6 0.8 I..0 

r/r, 

Fig. 13a. lonxza-txon Structure of Hydrogen and Helxum. 

The neutral fraction of hydrogen (solid lines) and 
helium (dashed lines) and (dash-dot) for uniform 
dustless nehulae as a function of normalized xadius 
-r/r^. The open circles denote the volume average of 

C^. 0.1. 
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Fig. 13b. Ionization Structure of Hydrogen and Helium. 

The neutral fraction of hydrogen and helium and Q 
as functions of t/Tj^ for uniform dusty nebulae with 
a constant dust cross section (O! := 0) , i=a 1.0 and 
Y = 0.1 for all cases. The symbols are the same as 
in Fig. 13a. 




Cg slowly but IS no longer negligible compared to unity. 

Therefore, neglecting ^2 commonly done, will result in an 

incorrect determination of the value of Yq'^^ from observations . 

For example,, for a given value of Y and observed value of 

(r^ IS the H Stromgren radius and is the He Stromgren radius), 

the required value of Yq underestimated when is neglected. 

In general, for small values of ^2'^^!^ ^2 negligible and is 

rarely exceeds 0.4. 

For nebulae with dust and with hydrogen and helium, simple analytic 
solutions to the coupled differential equations (B.14) are possible 
only if Ct IS negligible or zero, i.e. only if ct^(v) varies slowly 
so that a, ^ cr , In this case we have plotted on Figure 2-2. the 

A 

variation of Y with S (dashed lines). ITsing- these results, we 

J}. 

have calculated (1-x), (1-y) and' for uniform nebulae, e = 2 x 10 , 

and «« 1.0, These results are presented in Figure 13b. Comparison 
of Figures 13a and 13b show that for a given size r^ of the nebula, 
introducing dust increases the fractional ionization of both hydrogen 
and helium. Ihis is primarily due to the increased value of S^, and 
thus of S('t) throughout the nebula. 

TOien 0! 0, simple analytic solutions are not possible. However, 

Eq. (B.14) can readily be solved* numerically. On. Figure 12 we present 

A 

Y versus S for a few values of CC. In general we find that the 
shapes of the (1-x) and (1-y) curves are nearly independent of the 
details of the problem and are determined primarily by the- values of 

A 

Sq and Yq« Consequently, we do not present the xonization structures 
for a ^ 0. 


12k 



B.5 Observable Parameters 

Since we treat fluxes integrated over frequency, we lose most 
of the information on the variation of the spectrum of UV radiation 
throughout the nebula. Consequently, we cannot calculate the ioniza- 
tion structures of the trace elements with these semi-analytic models. 
The main result of these solutions is the ionization structure of 
hydrogen and helium. Hence we can calculate ^gas ^e have 

neglected doubly ionized helium so refers only to He^ ) 

The fraction of ionizing radiation absorbed ("destroyed”, cf. 
section 3.8) by gas is 


gas 


®0 


r. 

/ 


2 2 , 


dr(x + Yy)[X - Yy(l - n (r) 


0 


= S 


/I 

0 


(B.31) 


dt[x+ Yy(l - Ct)](x+ Yy)R ( t)G f t) /G .( t) 


For Eq. (b. 31) reduces to 




(B.32) 


where 


T 

■/■ 

0 


f(T) =/ dtR^(t)Gg(t)/G^(t) 


(B.33) 


and Tg dust optical depth at the helium Stromgren radius. 

The ratio of the intensity of helium to hydrogen recombination 
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lines (Eq. 3.36) is equal to 


(1+ Y)£(.2)/Sgf 


gas 


(b.34) 


As explained previously (Petrosian 1973^ 197^) ^ ^ given value of 

Y /Y the presence of dust with a,(v) constant (i.e. o: = O) 

'0 d 

increases the value of with increasing values of the total dust 

optical depth For negative values of Ct ^ ^ ^d 2 }^ fewer 

He ionizing photons compared with H ionizing photons are absorbed 
by the dust and is larger The reverse occurs for positive 

values of a. It should be noted the depends on Yq/^ 

insensitive to the value of Y. is also fairly independent of 

non-unif ormities or inhomogeneities in the dust and gas distribution 
as long as the dust-to-gas ratio is constant. However;, changes 

for dust and gas distributions that are not the same. 

Note that because we have used the approximations x = 1 for 


r s r. 


and y = 1 for r ^ the ratio 


S 1. 


126 



REFERENCES 


Aldrovandi, S. M. V. and Pequxgnot, D. , 1973, Astr, Ap, , 2 ^, 137, 

Auer, L. H. , 1967, Ap. J. (Letters), 1^, L53. 

Auer, L. H. , 1971, J. Quant. Spectrosc. and Rad, Trans,, 1^, 573. 

Auer, L. H, and Mihalas, D, , 1970, M.N.R.A.S,, 3^, 65, 

Auer, L.H, and Mihalas, D. , 1972, Ap, J, Suppl., No. 205, 193. 

Balick, B., 1975, Ap. J. , 201 , 705. 

Balxck, B. and Sneder, C. , 1976, Ap, J, , 208 , 336. 

Bohlin, R. C., Marxonni, P. A. and Stecher, T. P, , 1975, Ap. J. , 202 , 
415. 

Burgess, A, and Seaton, M. J. , 1960, M.N.R.A.S., 121 , 471. 

Chandrasekhar, S., 1960, Radiative Transfer (New York: Dover) p. 23. 

Dana, R. A. and Petrosian, V., 1976, Ap, J, , 208 , 354 (Paper II). 

Davison, B. , 1957, Neutron Transport Theory (Oxford* Clarendon Press). 

Dycke, H. M. and Simon, T. , 1977, Ap. J. , 211 , 421. 

Flower, D. R. , 1968, Planetary Nebulae (lAU Symposium No. 34), ed. 

D. E. Osterbrock and C.i R, O'Dell (Dordrecht; Reidel), p. 205, 

Greenspan, D, , 1970, Introduction to Numerical Analysis and Applications 
(Chicago* Markham), p. 53. 

Henrey, R. J. W. , 1970, Ap. J. , 1153, 

Hjellming, R. M. , 1966, Ap. J. , 143 , 420. 

Hummer, D. G. and Rybicki, G. B, , 1971, M.N.R A.S., 152 , 1. 

Hummer, D. G. and Seaton, M. J. , 1963, M.N.R.A.S., 125 , 437. 

Hummer, D. G. and Seaton, M. J. , 1964, M.N.R.A.S,, 127 . 217. 

Jameson, R, F. , Longmore, A. J. , McLinn, J. A. and Woolf, N. J. , 1974, 

Ap. J. , 190 . 353. 

Kirkpatrick, R.’ C. , 1972, Ap, J. , 176 . 381. 

Kirkpatrick, R, C. , 1977, private communication. 


127 



Kleinmann, D. E. , Gillett, F. C. and Wright, E. L. , 1976, Ap. J. , 208 , 

42. 

Knacke, R. F-. and Thomson, R, K. , 1973, Pub. A.S.P. , 8^, 341. 

Krugel, E. , 1975, Astr. Ap. , 38, 129, 

Lamarsh, J. R. , 1966, Nuclear Reactor Theory (Reading, Massachusetts 
Addison-Wesley) . 

Leung, C. M. , 1975, Ap. J, , 1 ^, 340. 

Mathis, J. S., 1970, Ap. J, , 1^, 263. 

Mathis, J. S., 1971, Ap. J. , 1^, 261, 

Mathis, J, S. , 1972, Ap. J. , 176 , 651. 

Mezger, P. G. , Smith, L. F. and Churchwell, E. , 1974, Astr. and Ap. , 32 , 

269, 

Natta, A. and Panagia, N. , 1976, Astr. Ap. , 191, 

Osterbrock, D. E. , 1974, Astrophysics of Gaseous Nebulae (San Francisco: 
Freeman) (AGN). 

Panagia, N. , 1973, A. J. , 78, 929. 

A 

Pequignot, D. , Aldrovandi, S. M. V. and Stasinska, G. , 1977, preprint. 

Petrosian, V. , 1973, in lAU Symposium 52 , Interstellar Dust and Related 
Topics , ed. J. M. Greenberg and H. C. van de Hulst (Dordrecht: 

Reidel), p, 445, 

Petrosian, V., 1974, in HI I Regions and the Galatic Center, Proc. 8th 

ESLAB Symposium , ed. A. F. M. Moorwood (Froscati: ESROSP-105) , p. 173, 

Petrosian, V., Silk, J. and Field, G. B. , 1972, Ap. J. (Letters), 177 , 

L69. 

Petrosian, V. and Dana, R, A., 1975, Ap. J. , 196 , 733 (Paper I), 

Rieke, G, H. and Low, F, J. , 1975, Ap, J. (Letters), 199 , L13. 

Righini, G. , Simon, M, and Joyce, R. R. , 1976, Ap. J. , 207 , 119. 

Rubin, R. H. , 1968, Ap. J. , 153 , 761, 

Sarazin, C. L. , 1977, Ap. J, , 211 , 772, 

Seaton, M. J. , 1960, Rept. Progr, Phys. , 313. 


128 



Simon, T. and Dycke, H. M. , 1975, M.N.R.A.S., 172 , 19P. 

Spitzer, J. , Jr,, 1968, Diffuse Matter in Space (New York* Interscience). 

Tarter, C. B. , 1967, Hi.D. Thesis, Cornell University. 

Telesco, C. M. , Harper, D. A. and Loewenstein, R. F. , 1976, Ap. J. 
(Letters), 203 , L53. 

Telesco, C. M. and Harper, D. A., 1977, Ap. J. , 211 , 475, 

Thronson, H. , 1977, Ph.D. Thesis, University of Chicago, 

Weisheit, J. C, , 1974, Ap. J. , 190 , 735. 

Weisheit, J. C. and Collins, L. A., 1976, Ap. J, , 210 , 299. 

Yorke, H. W. , 1977, paper presented at the Symposium on Recent Results in 
Infrared Astrophysics, January 1977, NASA-Ames Research Center. 


129 



